

DATl^PfiOCESSINp ALGOaXTRfPnJ RAL SCfiNNE * 

19 7 3 (EnvironmeD t al Research \ ns ± 1 of”' 
Hicfcigan)ff4* p HC $6.00 isci 05B 



NASA CR- /-? ^ 
ERIM 31650 -149-T 

N73-3034'71 


Onclas 

63/13 13181 


MULTISPECTRAL SCANNER DATA PROCESSING 
ALGORITHM DOCUMENTATION 


by 

M. Gordon and J. Erickson 
Infrared and Optics Division 


ENVIRONMENTAL 

RESEARCH INSTITUTE 
OF MICHIGAN 

FORMERLY WILLOW RUN LABORATORIES. 
THE UNIVERSITY O^' MICHIGAN 


Reproduced by 

NATIONAL TECHNICAL 
INFORMATION SERVICE 

U S Department of Commerce 
Springfield, VA. 22151 


prepared for 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Lyndon B. Johnson Space Center •> 

NAS 9-9784, Task B 2.14 




now* 


s»wi*o* 


IflWW® 


CUSVmIHS 

COLOR ILLUSTRATIONS 



NOTICES 


Sponsorship. The work reported herein was conducted by the Environ- 
mental Research Institute of Michigan* for the National Aeronautics and Space 
Administration, Lyndon B. Johnson Space Center, under contract NAS 9-9784. 
Dr. Andrew Potter (TF3) is Technical Monitor. Contracts and grants to the 
Institute for the support of sponsored research are administered through the 
Office of Contracts Administration. 

Disclaimers ^ This report was prepared as an account of Government- 
sponsored work. Neither the United States, nor the National Aeronautics 
and Space Administration (NASA), nor any person acting on behalf of NASA: 

(A) Makes any warranty or representation, expressed or implied with 
respect to the accuracy, completeness, or usefulness of the infor- 
mation contained in this report, or that the use of any information, 
apparatus, method, or process disclosed in this report may not in- 
fringe privately owned rights; or 

(B) Assumes any liabilities with respect to the use of, or for damages 
resulting from the use of any information, apparatus, method, or 
process disclosed in this report. 

As used above, "person acting on behalf of NASA" includes any employee or 
contractor of NASA, or employee of such contractor, to the extent that such 
employee or contractor of NASA or employee of such contractor prepares, 
disseminates, or provides access to any information pursuant to his employ- 
ment or contract with NASA, or his employment with such contractor. 

Availability Notice . Requests for copies of this report should be referred 
to: 

National Aeronautics and Space Administration 
Scientific and Technical Information Facility 
P. O. Box 33 

College Park, Maryland 20740 

Final Disposition . After this document has served its purpose, it may be 
destroyed. Please do not return it to the Environmental Research Institute of 
Michigan. 


* 

Environmental Research Institute of Michigan 
P.O. Box 618 

Ann Arbor, Michigan 48107 



NOTICE 

THIS DOCUMENT HAS BEEN REPRODUCED FROM THE 
BEST COPY FURNISHED US BY THE SPONSORING 
AGENCY. ALTHOUGH IT IS RECOGNIZED THAT CER- 
TAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RE- 
LEASED IN THE INTEREST OF MAKING AVAILABLE 
AS MUCH INFORMATION AS POSSIBLE. 



NASA CR- 
ERIM 31650-149-T 


Technical Report 


MULTISPECTRAL SCANNER DATA PROCESSING 
ALGORITHM DOCUMENTATION 


by 

M. Gordon and J. Erickson 
Infrared and Optics Division 


ENVIRONMENTAL 

RESEARCH INSTITUTE 
OF MICHIGAN 

FORMERLY WILLOW RUN LABORATORIES, 

THE UNIVERSITY OF MICHIGAN 

prepared for 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

July 1973 


NAS 9-9784, Task B 2.14 


Lyndon B. Johnson Space Center 
Earth Observations Division 
Houston, Texas 77058 


/ 



PRECEDING PAGE BLANK NOT FILMI© 



FORMERLY WILLOW RUN LABORATORIES. THE UNIVERSITY OF MICHIGAN 


FOREWORD 

This report describes pari of a comprehensive and continuing program of re- 
search concerned with advancing the state-of-the-art in remote sensing of the en- 
vironment from aircraft and satellites. The research is being carried out for the 
NASA's Lyndon B. Johnson Space Center, Houston, Texas, by the Environmental Re- 
search Institute of Michigan, formerly the Willow Run Laboratories of The University 
of Michigan. The basic objective of this multidisciplinary program is to develop re- 
mote sensing as a practical tool to provide the planner and decision-maker with 
extensive information quickly and economically. 

Timely information obtained by remote sensing can be important to such 
people as the farmer, the city planner, the conservationist, and others concerned 
with problems such as crop yield and disease, urban land studies and development, 
water pollution, and forest management. The scope of our program includes: (1) ex- 
tending the understanding of basic processes; (2) discovering new applications, de- 
veloping advanced remote-sensing systems, and improving automatic data processing 
to extract information in a useful form; and (3) assisting in data collection, pro- 
cessing, analysis, and ground -truth verification. 

The research described here was performed under NASA Contract NAS 9-9784, 
Task B 2.14 and covers the period from November 1, 1971 through January 31, 1973. 
Dr. Andrew Potter has been Project Manager. The program was directed by R. R. 
Legault, Associate Director of the Institute, and by J. D. Erickson, Principal Investi- 
gator and Head of the Multispectral Analysis Section. The Institute number for this 
report is 31650-149-T. 
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ABSTRACT 

The procedures followed in multispectral scanner data evaluation may be 
divided into two categories: (1) data handling and calibration, and (2) recogni- 
tion processing and subsequent evaluation of output. The report provides a 
description for the analyst of the algorithms employed in the current ERIM data 
processing scheme. Methods for suitable visual display of the results of this 
processing are also discussed. 
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MULTISPECTRAL SCANNER DATA PROCESSING 
ALGORITHM DOCUMENTATION 

l 

SUMMARY 

The classification and evaluation procedure of M-7 multispectral scanner data consists, in 
general, of two stages: pre-classification data handling and reformatting, and digital recogni- 
tion processing and evaluation. 

Pre -classification data handling and reformatting involves the following stage of proces- 

, sing: 

(1) Tape duplication, to make rough corrections for channel skew because of differences 
between recording and playback tape equipment 

(2) A-D conversion with "smoothing,” to convert analog data to digital data of a specified 
form while combining overlapping data and reducing the instrumental noise level of 
the data 

(3) Data misalignment correction 

(a) to evaluate the corrections necessary for channel skew, i.e., the misalignment of 
analog data on the physical tape because of improper alignment of the recording 
heads 

(b) to correct for scan line slew, i,e., the spatial misregistration of resolution ele- 
ments between scan lines 

(4) Clamping, scaling, and deskewing, to apply the channel deskewing corrections deter- 
mined by the program AUTOCAL, and subsequently to calibrate dynamically each scan 
line based on information present in each line 

The digital data tape is ready at this point for input into the classification programs. These 
classification procedures are: 

(1) Visual data display, to allow the user to select training set regions and the actual 
region to be classified 

(2) Signature extraction, to gather statistical information concerning the distribution of 
data values for a particular region in the target area in order to establish a decision 
criterion for the classification process 


1 



2 ™ 


FORMERLY WILLOW run LABORATORIES. THE UNIVERSITY OF MICHIGAN 


(3) Classification, to evaluate a set of resolution elements individually in terms of the 
signatures previously established, and associate each resolution element evaluated 
with one of the signatures, along with a measure of the likelihood of correct classifi- 
cation for each element 

(4) Recognition display, to provide a spatially registered visual display of the output of the 
classification procedure by means of printed characters, with features allowing the 
reliability of the recognition output data to be evaluated as well as the overall relia- 
bility of both the multispectral scanner data and collected statistical signature data 
used by the classification algorithm 

This set of procedures forms a sequence which, in most cases, generates an accurate repre- 
sentation of the distribution of substances and features in the target area of the scan. However, 
these procedures are not entirely automatic; at each stage the researcher must select various 
parameters necessary for proper execution of the algorithms defined in this report. 
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2 

INTRODUCTION 

Several basic algorithms are implemented in the ERIM multispectral data handling, pro- 
cessing, and evaluation system. This report outlines for the analyst the basic principles of the 
system without attempting a rigorous justification of the algorithms currently implemented. 

Although other types of multispectral scanner data (such as ERTS scanner data) are pro- 
cessed by ERIM, most of the data to be evaluated are obtained from the M-7 scanner [1|. Since 
data obtained from this instrument are in analog form, a set of manipulations must be performed 
to put them into a format acceptable to the equipment used in actual processing and evaluation. 
This process is known as the A-D (analog-digital) conversion process. A-D conversion is not 
necessary for out-of-house data, since in most instances this information is already in digital 
form acceptable to our bulk-processing digital computer (a C DC 1604- b). 

Section 3 deals with M-7 scanner data collection procedures, subsequent A-D techniques 
for reformatting the data, and corrections applied in the A-D process to compensate for pos- 
sible instrumental inconsistencies. Section 4 discusses the methods by which additional cor- 
rections, calculated by techniques discussed in Section 3.3, are applied to the data, as well as 
factors inherent to the data itself. After these data preparation processes are completed, vis- 
ual displays of the data are generated for use in evaluating the raw data and to aid the user in 
evaluating the best method available for futher data processing. These display methods are 
discussed in Section 5. Sections 6 and 7 deal with the methods by which the data is evaluated 
and classified, that is, the methods by which data points are grouped together and associated 
with a particular class which has been defined by a pre -determined distribution called a sig- 
nature. Finally, Section 8 deals with the methods by which the results of classification or rec- 
ognition processing are displayed. 
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3 

ANALOG-DIGITAL DATA CONVERSION 

3.1. M-7 SCANNER AND DATA COLLECTION 

The M-7 scanner, housed in a C-47 aircraft along with analog recording equipment for bulk 
data storage, is a device which senses the radiation from an area and, by means of a system of 
optics, divides this radiation into a number of spectral regions (channels) covering the visible 
through infrared regions of the spectrum. Figure 1 outlines the general configuration of the 
M-7 scanner system and peripheral equipment. 

A scan mirror, internal to the scanner housing, is rotated in a plane essentially perpendi- 
cular to the direction of flight (see Fig. lc), such that the sensors within the scanner are se- 
quentially viewing (1) the scene below, (2) one of several internal calibration sources, or (3) the 
internal scanner housing, which provides a measure of the dark level reading of the sensors. 
Each rotation of the scan mirror generates a set of data called a scan line and effectively con- 
tains a measure of the radiation emanating from a swath of target area, such that the largest 
dimension of this rectangular swath is perpendicular to the direction of flight. This informa- 
tion from each of the spectral channels, along with additional synchronization information, is re- 
corded continuously on the analog tape. After a desired set of successive scans (referred to 
as a run) has been recorded on an analog tape, the tape is brought, to the ERIM data processing 
facility for digitization and subsequent processing. 

3.2. ANALOG DATA AND TAPE DUPLICATION 

The M-7 multispectral scanner data is recorded on analog tape and must go through A-D 
conversion to produce a digital equivalent compatible with the CDC 1604 -B computer. 

Because of the large quantity of data recorded on the analog tape, only certain areas of the 
total scan line are usually digitized. A CRT x-y display is used for a C-scan presentation of 
one channel at a time in the video region of the analog tape (the region wherethe sensor was actu- 
ally viewing the scene and notsome internal portion of the scanner). This C-scanpresentation 
essentially provides a visual display of the radiation recorded by the scanner, making gross 
target features apparent. As this C-scan is being presented, a line count is also displayed so 
that scan line numbers for certain regions of interest can be identified for use in the digitizing 
procedure. 

As the scan mirror rotates, a separate channel of information containing two synchroniza- 
tion pulses per scan line is recorded on the analog tape (Fig. Id). One of these pulses, the 
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FIGURE 1. M-7 SCANNER CONFIGURATION 
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FIGURE 1. M-7 SCANNER CONFIGURATION (Continued) 
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(c) Geometry of Airborne Scanning 
FIGURE 1, M-7 SCANNER CONFIGURATION (Continued) 
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marker sync pulse, denotes the time when the rotating scan mirror passes a fixed point rela- 
tive to the scanner housing. The other sync pulse, called the roll stabilized pulse, indicates the 
time when the scan mirror passes a fixed point relative to the nadir. From the difference in 
the occurrence time of these two pulses, it is possible to measure the amount of aircraft roll from 
scan line to scan line. These differences may be digitized along with the multispectral data 
with special equipment that generates a separate channel called the time channel. 

The tape recorder heads used to record the analog data in the aircraft generally have a 
different orientation relative to the direction of tape movement than the tape recorder used to 
play back the analog tape to the A-D converter. Because of this, a duplicate of the original 
analog data tape is made that introduces a time lag, usually different for each track, by means 
of variable delay transmission lines. This delay process allows more exact registration between 
data channels. The skew between any two channels is usually reduced to less than 2 ju sec by 
this procedure. In theory, it is possible to correct for skew between channels by measurement 
of the differences in head alignment between the recording and playback equipment. In practice, 
however, the amount of time required for these measurements is prohibitive, and therefore the 
amount of delay to be applied to each channel during tape duplication is only approximated. 

Any further channel skew is calculated and corrected after the A-D conversion process is com- 
pleted by the program AUTOCAL. discussed in Section 3.4. The corrected duplicate copy of the 
analog data tape is then used as input to the A-D processor for conversion to digital form. 

Several regions along each scan must be digitized, because they contain information nec- 
essary to calibrate and cross-correlate the absolute data values in the video region, once these 
are in digital form. Built into the A-D converter is a resolution rate generator, which divides 
the region between two successive marker sync pulses into a specified number of regions each 
of which is digitized separately. However, it is not necessary to digitize the whole region be- 
tween two sync pulses, since much of this is merely the scanner viewing the opaque inner wall 
of the Instrument. The A-D converter may be gated manually such that only certain regions are 
digitized. The regions generally digitized are: 

(1) video, that region where the sensor views the scene below the aircraft 

(2) dark level, a portion of that region where the sensor views the dark interior of the 
scanner housing 

(3) sun sensor, that region where the sensor views the radiation passing through a diffuse 
opal glass plate on the top of the aircraft 

(4) calibration lamp, that region where the sensor views a radiance -transfer standard 

(5) cold plate, that region where the sensor views a source of cool, known temperature 

(6) hot plate, that region where the sensor views a source of known temperature, hotter 
than that of the cold plate 
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Figure 2 shows the time correlation between the analog data in a particular channel and 
the sync pulse channel, the digitizing gates, and the resolution rate generator pulses. 

3.3. ANALOG-TO-DIGITAL CONVERSION WITH SMOOTHING 

The actual digitizing procedure performed by the A-D converter is fairly simple: 

(1) A resolution element which has been selected to be digitized by means of the gate 

settings is read by the analog tape recorder; the data value of each channel is strobed 

into a 48 -bit buffer and shifted right the appropriate number of places as designated 

by a data resolution switch containing the number of significant bits to be used per 

data value. The data resolution switch is ordinarily set manually to 9 bits per data 

_q 

element, resulting in an accuracy of 2 over the analog voltage range. 

(2) The completion of a 48 -bit word causes the computer to transfer that word to storage 
in sequential, increasing addresses in core. 

(3) The CDC 1604 computer is signalled when the final resolution element for a particular 
scan line has been digitized and transferred to the CDC 1604 computer; a flag is sent 
to the 1604 that the block of data is complete. The computer then writes the digitized 
information onto magnetic tape, one scan line per digital tape record. (Appendix I 
describes the format of the digital tape.) 

Flight altitudes of 2000 ft and above ensure some overlap in successive scan lines, i.e., 
several successive scan lines contain data from portions of the same region on the ground. 

This results in a quantity of analog data much larger than is necessary for reasonable data 
processing and evaluation. Thus, as the analog-digital conversion is taking place, an on-line 
computer program, designated A2F3, is being executed on the digital computer. This program 
performs a function called smoothing, or filtering, which has a twofold purpose: 

(1) The reduction of the number of scan lines by a factor between 3:1 and 16:1. 

(2) The filtering (or reduction) of instrumental noise present in the recording process. 
This is accomplished by simple averaging of each data value over NSMOOTH lines, where 

NSMOOTH is the number of lines to be combined to form a single digitized line. If the dis- 
tribution of instrumental noise associated with each scan line is assumed to be Gaussian, 
averaging cancels out much of this background. 

The actual computational procedure uses the following areas within the computer: 

(1) Input buffer of N WORDS locations where 
N WORDS - (NSS * NCR AN + 4)/5 

NSS = the number of resolution elements per scan line 
NCHAN = the number of channels per resolution element 

(2) Scratch buffer consisting of two memory locations, used for intermediate, partial un- 
packing 10 
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FIGURE 2, SCAN LINE WITH ENABLE GATES AND GATED SAMPLE PULSES 
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(3) Two averaging buffers consisting of NWORDS memory locations, each of which is 
used to accumulate the data value averages 

(4) Output buffer of NWORDS memory locations, used for final repacking and output to 
magnetic tape 

* A memory location on the CDC 1604-B computer is composed of 48 -bits; as scan lines are 
converted from analog to digital form they are directed to the computer memory and stored, 

5 data values of 9 bits per word, with 3 bits unused (Fig. 3a), These data values are in 2's 
complement notation. Complement notation is described in Appendix I. The computational 
procedure used is: 

(1) The first of NSMOOTH scan lines is read into the input buffer. Locations 1 through 
NWORDS are successively subjected to the following procedure: 

(a) Each of the 5 data values in location N (where 1<N<NW0RDS) is biased upward by 
complementing its sign bit, 

(b) Each value is partially unpacked into one of two locations called partial unpacking 
buffers, the even samples into one location, the odd into the other. (See Fig. 3b.) 

(c) If NSMOOTH > 8, the data values are each divided by 2 and the remainder dis- 
carded, with an accompanying loss of significance. 

(d) The contents of these two locations are then stored in one of two (even or odd) 
averaging buffers. 

(2) Scan lines 2 through NSMOOTH are then successively read in, with steps la through 
Id applied to each location of each scan line; however, step Id is now modified 
such that the partial unpacking buffers are not stored in the averaging buffers, but are 
arithmetically added to them. 

(3) After scan line NSMOOTH has been processed, the averages are truncated to 9 bits 
per data value and multiplied by 1/NSMOOTH, where this multiplier has been pre- 
viously truncated to 9 bits. The product is then truncated to 9 bits. 

(4) The odd and even data values are repacked into the output buffer in a slightly different 
format (Fig. 3c). the bias removed by recomplementing the leftmost bit of each data 
value, and the data values converted to l's complement by adding +1 to all negative 
data values^ 

(5) Finally, the output buffer is written in ADTEST2 format (see Appendix n) on the output 
magnetic tape. 
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This procedure results in no loss of significance for smooths of 4- and 8 -to -1; however, 
smooths of 9- through 16-to-l initially lose 1 bit of significance, and smooths other than 4-, 

8-, or 16-to-l lose 1 bit of significance at the end of averaging. 

3.4. AUTOCAL— EVALUATION OF CHANNEL SKEW AND SCAN LINE SLEW 

After the A-D conversion process is completed, the resultant digital tape is used as input 
to a CDC 1604-B computer program, called AUTOCAL, which has a three-fold purpose: 

(1) Calculation of an average scan line, corrected for slew effects (spatial misregistration) 
between scan lines 

(2) From the calculated average scan line, a determination of the amount of skew present 
between channels, based on the position of the sun sensor peak in each channel (see 
Section 3.2) 

(3) Calculation of a dark level value for each channel of data 
THE ALGORITHM SEQUENCE 

The user identifies the channel (ICHAN) to be used as a reference in the calculation of 
skews for the other channels, and the boundaries Cl and C2 defining the correlation region 
used in both the deslewing and deskewing procedures. The program individually processes 
successive scan lines by editing noise spikes from the data, correcting for slew, and accumulat- 
ing the sum line. 

SPIKE NOISE EDITING 

Data differences — i.e«, the differences between the data values in each channel for successive 
resolution elements along the scan line — are computed. A noise spike is characterized by two 
very large differences of opposite sign enclosed by two ordinary differences. If the minimum 
of the two large differences is greater than 5 times the maximum of its two neighboring dif- 
ferences and also greater than a high estimate of the expected difference, DEL(J), a noise spike 
is identified, and the data value then replaced by the average of the preceding and succeeding 
values in the same channel. DEL(J) is obtained from the first line by ranking the absolute dif- 
ferences, and then multiplying the 80-th percentile difference by 2. DEL(J} is included to cover 
the case in which both the neighboring differences are very small or even zero. A steeply rising 
oi falling signal is not identified as a noise spike if the neighboring pulses are much the same 
size and if the middle differences are of the same sign. 

SLEW CORRECTION 

Assuming that the scan lines consist of N resolution elements and NCHAN channels per 
resolution element, data differences (as defined previously) in channel ICHAN from the first 
data line between points Cl and C2 are computed and stored in a vector Y. Channel ICHAN is 
defined as the reference channel; Cl and C2 have been previously defined; and Cl < C2 < N. 
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Thus Y is a vector containing elements through Y^ where C - C2 - Cl + 1. Then a set of 
sums, S, is formed such that 

j+C+1 

S i " T. P k.ICHAN Y k-j+l 
M 

where i = 1, N - C + 1 and where D k ICHAN is defined as the data difference between the k-th 
and (k + l)-st resolution element in channel ICHAN for the second scan line. Since maximum 
correlation occurs when this sum is a maximum, the data line is shifted in such a way as to 
align this region of maximum correlation with the region from Cl to C2 on the first scan line. 

(If there were no slew, the region of maximum correlation would occur in the Cl to C2 region 
for each scan line.) 

After all the scan lines have been processed, the program computes the average line in 
each channel by dividing each element of the sum line by the number of lines processed. This 
average line is used for the calculation of the skew of each channel. The standards are points 
Cl to C2 on the ICHAN average line. Each of the other lines is tested to find the region in 
which the length C2 - Cl + 1 correlates most closely with the standard. The method of correla- 
tion is, as before, the sum of products of differences, but instead of locating a local maximum, 
all possible sums with lags of -10 to +10 are tested and an absolute maximum obtained. The 
result of the correlation is output as a set of integers showing the number of points by which 
each channel precedes (negative) or follows (positive) the reference channel. Also, fractional 
slews for each channel are estimated by passing a parabola through the biggest sum and the 
two neighboring sums. The location of the base of the parabola determines the fractional slew, 
which differs from the integer slew by no more than 0.5, a figure achieved when two adjoining 
sums are equal and maximal. 

Then a sum line, A. is cumulated from the first data line and the shifted data line, such 

that 


N NCHAN 

\rE L 

i.i j.i 




where is the first scan line and is the second scan line. 

Processing of this kind proceeds with successive lines (all lines being slew “Corrected with 
reference to the first scan line). Each new shifted scan line is cumulated with the sum line A, 
such that 
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N NCHAN 

\r£ £ “ij'Vi 1 

i-1 j=l 


where V represents the m-th data line alter shifting for slew correction, 
m 

Finally, a dark level for each channel is computed from the average line. This dark level 
for channel j is defined as 

DARK. = MAX(S.) 


i+7 

where S. = A. . for i = 1 to N - 7 
k=i 


Each value of the average line is subtracted from this DARK array, generating a corrected 
average line C, such that 

0. . - DARK. - A. . where i = 1 to N 
M J i,] 

and j = 1 to NCHAN 


The output from AUTOC AL provides information which may be used to provide further chan- 
nel registration corrections and better boundaries for calibration regions, as inputs to the next 
stage of digital processing: clamping, scaling, and deskew (namely, the program CSD). 


4 

CLAMPING, SCALING, DESKEWING 

Any or all of three types of corrections to the data contained in each scan line may be 
needed after an analog tape has been digitized: 

(1) dark level correction 

(2) multiplicative scaling based on one of the calibration sources 

(3) deskewing of the channels to correct spatial misalignment 

Dark level correction and calibration scaling both employ data values contained in each 
scan line to calculate the correction factors. The deskew correction factors obtained from the 
AUTOCAL program are now used as input to this correction program, called CSD. 

The user must supply the following correction data to the program CSD for the generation 
of a new data tape: 

(1) The line numbers of the scans to be corrected and the point numbers which represent 
the video portion of these scans 
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(2) N skew values where N Is the number of data channels on the tape, and the reference 
channel is assumed to have a skew of zero 

(3) For clamping (dark level correction), either 

(a) a clamp region — i.e., line and point numbers of a dark level region, from which a 
dark level value for each channel is calculated 

(b) actual clamp levels (dark levels) to be used 

(4) For scaling 

(a) the point numbers of the region for the calibration signal of each line 

(b) the sampling region line numbers if time channel sampling (Section 3.2) is to be 
used 

(c) a smoothing constant representing the number of lines over which the peak in the 
calibration region is filtered by means of an exponential-type function 

(d) multiplicative scale factors (optional) 

4.1. CLAMPING 

The user either specifies a clamping region to set the dark levels of each individual scan 
line, or supplies a constant dark level value for each channel. 

If the dynamic dark level correction method is selected, the points specified by the user 
within the clamping region are averaged channel by channel, and an array DARK(I), 1=1, 
NCHAN is generated for each line. If dynamic dark level correction is not used, DARK(I), 

1 = 1, NCHAN remains constant throughout the run. 

The correction is applied such that for DATA(N, J) where N = 1, NSS and J = 1 NCHAN is 
modified to read DATA(N, J) = DATA(N, J) - DARK(J). 

4.2. SCALING 

Three scaling options are available: 

(1) dynamic scaling with a time channel criterion 

(2) dynamic scaling without a time channel 

(3) multiplicative constant scaling (user -specified) 

The use of any of these options alone is possible, as is a combination of either 1 and 3 or 

2 and 3. 
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If a time channel (Section 3.2) is available, the procedure is as follows: 

(1) A region is specified by the user, and the range of values for that region is calculated. 

(2) Each line of the entire scene to be processed is then tested to see if the value of the 
time channel for the particular line is within the range of time channel values for 
the region specified by the user in (1) above. If the value for the particular scan 
line is within this range, dynamic scale factors are calculated as for scaling option 
(2) above and applied to the scan line. If the time channel value for the scan line is 
outside the acceptable range (as defined by the user -specified region), the scale fac- 
tors calculated for the previous line are applied to the current scan line. 

Dynamic scale factors are calculated as follows: 

(1) A value is calculated for each resolution element of the scaling region as specified by 
the user; this value is the sum of all the channel values for the resolution element. 

(2) The maximum of these values is found, and data values in all channels at this peak 
element are filtered, by means of an exponential function to avoid noise spike problems, 
with the previous LSMOOTH lines where LSMOOTH is the number of user-specified 
lines for the smoothing process (default - 60), 

(3) The smoothed values for each channel are then used to divide the data in the video 
region. 

4.3. DESKEWING 

For NSS resolution elements per scan line in the video region and NCHAN channels per 
resolution element, the user must supply NCHAN values into an array ISKEW(I) where 1 = 1, 

10 (1— I — 10 in increments of 1). The variable MINSK EW is set equal to the minimum value 
of all of the NCHAN skew values in the ISKEW array. 

Then the ISKEW array is modified such that MINSKEW is subtracted from all of the NCHAN 
values of the ISKEW array. 

when DATA(N, I) (the value of the J-th channel for the N-th resolution element where J = 1, 
NCHAN and N - 1, NSS) is stored for scaling, DATA(N, J) is not returned, but replaced by the 
value DATA [N + ISKEW(J), J] . 
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DATA DISPLAY 

Displays of the multispectral data are produced by means of the CDC 1604-B computer 
program GRAY2. 

The user specifies a method by which the data values for an array called LEVEL are estab- 
lished. The LEVEL array contains data values which encompass the entire range of possible 
values for a given set of data (i.e., 0-511 for positive data, and -255 - 255 for bipolar data). 
These levels increase in value with increasing subscript, and a set of 2 characters is associated 
with a given level. 

The following methods are used to perform level settings: 

(1) Manual method: The user feeds levels directly to the program. The program per- 

forms elementary tests to ensure that the level values are monotonically increasing 
and encompass the entire range of possible data values. 

(2) Automatic level set by random sampling: This method randomly selects points from 
a user-specified region (usually a subset of the region to be displayed) using a pro- 
grammed random number generation routine. Then NLEVEL display levels (NLEVEL 
being a user-defined variable specifying the number of display levels) are set such 
that the distribution of the randomly- selected data values is uniform in each range 
delineated by the LEVEL values. 

(3) MINMAX method: A minimum and a maximum range value, MIN and MAX respectively, 
are user-specified, along with NLEVEL; levels are generated such that the difference 
between values of the LEVEL array is DIFF = (MAX - MIN)/NLEVEL. 

When the levels have been determined, association of each resolution element with one of 
the NLEVEL levels takes place in one of three modes. 

Simple classification : The data value in the channel to be displayed for a given resolution 
element is compared to the i-th level where i = 1, NLEVEL. When the data value for the resolu- 
tion element is greater than LEVEL(I - 1) but less than or equal to LEVEL(I), the resolution 
element is assigned the set of symbols associated with the f-th level for printing. 

Two-channel interval criterion classification : The level to which a particular resolution 
element is classified is selected In the same manner as for the simple classification method; 
however, the value of the resolution element in a second channel (a control channel) is then 
tested. If the value in the control channel is less than a user-supplied maximum value, the 
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resolution is assigned the I-th set of symbols for printing; however, if the value in the control 
channel is larger than the maximum allowed, the resolution is assigned a blank character for 
printing. 

Two-channel level-equality criterion classification : The level to which a particular resolu- 
tion element is classified is selected in the same manner as for the simple classification method; 
however, the value of the resolution element in a second channel (a control channel) is compared 
to the values of the array LEVEL. If the control channel value matches one of the LEVEL values, 
the resolution element is assigned the I-th set of symbols (where the element was classified as 
the I-th level initially); otherwise, the resolution is assigned a blank set of symbols for printing. 

Each point of the scan line is processed in this manner (that is, associated with a particular 
level and thus with the particular set of print symbols characteristic of that level); then a com- 
puter printout is produced, such that one line of printed symbols bears a one-to-one correspon- 
dence to the data values for a given scan line. Figure 4 exemplifies the format of this visual 
multispectral data display. A discussion of recognition color mapping is contained in Section 7. 
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FIGURE 4. EXAMPLE RECOGNITION DISPLAY 
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6 

SIGNATURE EXTRACTION 

The end goal of multispectral data processing is determination of the object class of a 
given resolution element in the target area. This is accomplished by defining areas, called 
training sets, which are assumed to be homogeneous and representative of a particular object 
class. The data values in a training set region are combined such that the set of mean and co- 
variance matrices, along with secondary statistical information, define what is called the sig- 
nature of a particular training set. Once a representative number of training sets and their 
associated signatures have been obtained, classification procedures (the assignment of a given 
resolution element to a particular class, based on comparison of the resolution element’s data 
values with the set of signatures) are instituted on the resolution elements of a particular sub- 
set of the data, called the scene, chosen by the user. 

This section will now describe the algorithms implemented for the extraction of signatures 
by the CDC 1604-B program SIG1. 

In the calculation of a given signature, the algorithm implemented in SIG1 is as follows: 

(a) All training sets must be assumed rectangular, and the user must supply scan line and 
point boundaries for each training set. 

(b) The mean value in each channel, a covariance matrix, the standard deviation from the 
mean in each channel, and a correlation matrix, are accumulated by processing the 
points within the training set individually. 

Given a two-dimensional DATA array, where DATA (I, J) refers to the value of the J-th 
channel of the I-th point to be processed, the processing of n data points will result in: 
n 

MEAN(J) = J ) 

1=1 

For the covariance matrix entries for channels J and K (where J and K are less than 
or equal to NCHAN)jthe total number of channels used to calculate the signature is 


COV(J, K) = 



L>1 


DATA(I, J) * DATA(I, K) 


- MEAN(J) * MEAN(K) 
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the standard deviations 


o T = 7COV(J, J) 

J 

and the correlation matrix 


COR(J, K) = 

J K 

(c) The values calculated for each training set are based on n data points being used. The 
program edits data points within the training set region so that extreme points are not 
used in the signature calculation. 

The editing procedure begins with the establishment of editing criteria. Each component 
of the data point is compared to the edit values. If any component is outside the acceptable 
range then the whole data point is rejected. 

The establishment of editing criteria requires an estimate of the parameters of the distri- 
bution being used so that extreme points or points not from that distribution may be excluded. 

A criterion based on a sample median and sample quartiles is preferable to one based on a 
sample mean and sample standard deviation because extreme points have much less effect on 
the median and quartiles. With the assumption of a form for the underlying population distri- 
bution and a probability threshold for exclusion, one can set upper and lower bounds using the 
median and quartiles. An underlying multivariate normal population is assumed, so the median 
as an estimate of the mean and the quartiles can be used to estimate the standard deviation. 

To use the normal (z) distribution, it is necessary to assume that the median is the popula- 
tion mean value and that the average quartile deviation from the median represents a value 
z_ - 0.6745 (see Fig. 5). Thus, 


-l/2(X 3/4 - X l/4 n x ( 


Q 


_ 1IQ 

a 


so 


a = 


X 


Q 


0.6745 


Next, a probability threshold criterion for excluding points is picked, based on the above 
assumptions. For SIG1, only one point in a thousand is to be rejected if the data meet the as- 
sumptions. This is equivalent to saying that the point should be rejected if 
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X - Median 
a 


3.2905 (See Fig. 6) 


The SIG1 edit routine considers the individual channel values to be independent normal 
samples, computes bounds separately for each channel, and rejects an observation point if any 
one of its channel values exceeds its bounds. Thus, the probability of any one channel exceeding 
its bounds must be (l/NC) times the criterion value (0.001), where NC is the number of channels. 
Thus, if NC - 10, so P(UI > z^) = 0.0001, then the criterion is % T = 3.89. Since the threshold 
X values, X^, are given by; 


X T (NC) = Median ± z t (NC)ct 


- Median ± 


z T (NCr 

L Z Q . 


X. 


z t (NC) 

= Median * 0.6745 X Q 
We have, for NC = 10, 

X T (NC) = Median ± 5,767 X Q 

Note that there is no dependence on the number of samples used, except indirectly in that 
the median and quartiles will more accurately estimate the true population parameters as more 
samples are used. 

The algorithm employed to determine the editing criteria is as follows. 

The first N data points in the area being processed are used to form the sample population. 
The value of N depends on the number of channels being used; N is given by the formula: 

N = 2400/NC 

also 

0 < N < 600 


Computer storage limitations account for the limit of an upper bound on N. 

Then, for each channel, the median and quartile deviations are found. The editing bounds 
are calculated as: 
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FIGURE 5. 


NORMAL DISTRIBUTION. 


Quartile deviation. 



FIGURE 6. REJECTION CRITERIA 
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HI(I) = MEDIAN + QDEV * EDFACTOR(NC) 


LO(I) - MEDIAN - QDEV * EDFACTOR(NC) 

rz (Ncr 

where EDFACTOR(NC) is the p 067 45 t errn from before. 


The values for EDF ACTOR for channels 1 - 24 are: 

CHAN 12345678 


EDFCTR 

4.88 

5.16 

5.32 

5.43 

5.51 

5.58 

5.64 

5.69 

CHAN 

EDFCTR 

9 

5.73 

10 

5.77 

11 

5.80 

12 

5.83 

13 

5.86 

14 

5.89 

15 

5.91 

16 

5.94 

CHAN 

EDFCTR 

17 

5.96 

18 

5.98 

19 

5.99 

20 

6.01 

21 

6.03 

22 

6.05 

23 

6.06 

24 

6.08 


Where LO(J) and HI(J) are the lower and upper editing limits, respectively, in the J-th 
channel, the individual points are not deleted if 

LO(J) £ DATA(I, J) < HI(J) 

for all channels, where DATA(I, J) is the value of the J-th channel for the I-th resolution ele- 
ment within the training set. 
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7 

RECOGNITION PROCESSING 

Once signatures have been extracted from a scene of multispectral data, it is then possible 
to attempt to associate a given resolution element with one of the known signatures. This pro- 
cess is called classification of multispectral scanner data, and is performed by the CDC 1604-B 
program EXPMAP. 

With M signatures availables for a given scene of data to be classified, a value Q is com- 
puted for each of the M signatures for each resolution element to be processed, where 

Q. =(X - Mi ) T R 1 “ 1 (X - M t ) 


where X = the vector containing the data values in all the channels used 
jll = the mean vector for the i-th signature 
R. - the covariance matrix for the i-th signature 

The resolution element in question is said to be classified as the j-th signature class, when 

Z , where Z. = Q + log I R, I, is greater than the Z computed for any of the other signatures. 

J J J J K 

Initially, before any classification processing takes place, the user must input an array S 

such that the k-th signature is associated with the element S^, where 0 < < 511 for all k. As 

the i-th resolution element is classified as belonging to signature k, a two-channel output scan 

line is being constructed such that the first channel contains the value and the second channel 

the value EXP., where 
r 

EXP. = Q.. * 5.12 + 0,5 when Q.. £ 99,6 

and is the Q value computed above for the i-th resolution element and the j -th signature class. 
If Q.j > 99.6, EXP. = 511. 

Thus the second channel can be used as a criterion denoting how well the classification 
procedure worked for each resolution element. 
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RECOGNITION MAPPING 

After the classification or recognition procedure has been completed, display of the results 
in a useful manner is necessary so that signature selection and classification techniques can be 
evaluated. Since the output from the EXPMAP classification program is another digital tape of 
the standard format (described in Appendix II), the display of the data can be accomplished by 
the program GRAY2 (Section 5). 

The two-channel modes described in Section 5 are used primarily for displaying recogni- 
tion output tapes. The two-channel interval criterion mode is used to display only those points 
whose classification was reasonably certain, based on the value stored in the second channel by 
EXPMAP (Section 7). It is possible with this feature to generate displays with different cut-off 
values for this second channel, and therefore with different error bounds. 

The two-channel level-equality criterion mode allows the user to evaluate the probability 
of correct classification over the scene for a given signature class. Consequently, the user 
might find that at the edges of a particular region the probability of correct classification was 
poorer than at the center of the region. 

The visual display of recognition output is extremely helpful in evaluating the recognition 
process, since it allows a large amount of data display in a relatively small amount of output 
in a form easily scanned by the user and in which gross features become more prominent than 
individual resolution elements. 

Figure 4 is an example of a color -coded recognition output visual display. 
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Appendix I 

COMPLEMENT NOTATION 


One's Complement 

The one’s complement of a binary number is defined as that number which results by re- 
placing each digit 1 in the original number with the digit zero and replacing each zero with 
a 1. For example, the one’s complement of the binary number 101 is the binary number 010. 

Note that the sum of a binary number and its one’s complement is the binary number composed 
of all ones (e.g., 101 + 010 = 111). 

Two's Complement 

The two's complement of a number is defined as the one’s complement increased by 1. 

Note that if an n -digit binary number is added to an m -digit binary number (where n £ m), any 
carry to the (n + l)st position is ignored. Thus a binary number added to its two’s complement 
is a number consisting of all zeroes (101 + 011 - 000). 
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Appendix II 
DIGITAL TAPE FORMAT 

One set of scan lines as recorded on analog tape results in a set of digital records, termi- 
nated by an end -of -file, as follows: 

(1) The first record, 2592 bits (54 CDC 1604 -B words) in length, is called the title record 
and contains information pertinent to the data set which will follow. The content of 
this title record is: 


1 

word 

Fortran code word (1777777777777777 8 ) 

12 

words 

BCD title of data set 

1 

word 

number of resolution elements per scan line (NSS) 

1 

word 

beginning angle from nadir (BANG) 

1 

word 

angular increment (DANG) 

1 

word 

number of channels (NCHAN) 

1 

word 

conversion factor (CONV) 

1 

word 

data positive flag (IPOS) 

1 

word 

data packed flag (I PACK) 

24 

words 

descriptive information in BCD 

10 

words 

unused 


(2) The rest are data records, containing one scan line of data per record. The number 
of CDC 1604-B words per record is 

NWORDS = p NSS * NCHAN) - lj + 1 

Each word contains five 9 -bit data values (Fig. 3), such that the data values representing 
channels 1 through NCHAN for the i-th resolution element are contiguous. 

(3) The data set is terminated by a standard end-of-file mark. 
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Appendix III 
PROGRAM LISTINGS 


PROGRAM S I G 1 

C NEEDS PROCESS# UNPACK3, POSDEF 

C VERSION 1.0 ( PROGRAMMERS W.L. BROWN 4 0.2UK# 12/711 

***W AR NI NG »* » SURR. PCSDEF USES 2»NC*»2 4 NC 4 3 SPACES IN ERASABLE 

•••FOR THIS REASON, THE 'DATA 1 ARRAY SHOULD BE KEPT FIRST IN COMMON 


'•* ' S I G1 ' IS A SPECIAL VERSION OF 'IMPROVE* WHICH ENABLES ONE TO CAL- 
CULATE STATISTICS FOR SEQUENTIAL# CONTIGUOUS TRAINING AREAS OF 
EQUAL SUE AS WELL AS 'NORMAL* TYPES OF TRAINING SETS 
'SIGl' IS A GENERAL SIGNATURE-CALCULATING PROGRAM WHICH REPLACES 
•CRECCAST' ... A FORM OF HISTOGRAM IS INCLUDED IN THE OUTPUT 
DOES NOT ALLOW TRANSFORMATIONS 


PACKED DATA TAPE ON 3, BCD OUTPUT ON A 

PUNCHED OUTPUT CAN BE WRITTEN ON ANY TAPE DRIVE OR ON 'O' TO DISCARD 
BCD OUTPUT OF TAPE CHANNEL [CP ON 5 ( ICP=0 TO OMIT OUTPUT ON 51 

INPUT MEDIUM IS T, F, 1, 2, 6, 7 OR 8 

•INPUT MEDIUM* IS REQUESTED PRIOR TO READING IN PROGRAM CONSTANTS 
AND THEN AGAIN PRIOR TO READING IN LINE# POINT DESIGNATIONS... THIS 
ENABLES ONE TC USE THE SAME CONTROL CARDS AS FOR * HIST1 • AFTER 
SECOND ASSIGNMENT OF INPUT MEDIUM IS MADE 

up to 13 channels# 6000 samples/scanline 

SENSE SWITCH 1 TO TURN OFF EOITING 

SENSE SWITCH 2 IF RUN CONTINUES OVER 2 OR MORE FILES 
SENSE SWITCH 3 FOR OPTION TO CHANGE TITLES 


FOR CHANNEL NOS. INSTEAD OF WAVELENGTH BANDS IN OUTPUT, 
•SPECTROMETER CHANNELS® ' WITH 0 OR WITH '15* 


ITROMETER IDENTIFICATIONS ARE 

• ♦ 

• 

SPECTROMETER CHANNEL 

BANDPASS 

1 

.43 


.47 

2 

.47 

- 

.49 

3 

.49 

- 

.51 

4 

.51 

- 

.53 

& 

.53 

- 

.56 

6 

.56 

- 

.59 

7 

.59 

- 

. 63 

8 

.63 

- 

. 67 

9 

.70 

- 

.90 

10 

1.00 

- 

1.40 

11 

1,50 

- 

1. eo 

12 

2.00 

- 

2.60 

13 

9.30 

- 

11.7 

15 

NO 

LABEL 


ANSWER 


C (BANDS ARE FOR M-7 SCANNER. CF. E. WORK 7/6/71 MEMO) 
C NC = NO. OF CHANNELS IN SUBSET. (OPTIONAL) 

C SUBSETS ARE READ BY (2013) 

C ICP = TAPE CHANNEL TO BE PLOTTED ON TAPE UNIT 5 
C ICP = 0 BYPASSES WRITING ON TAPE 5 


C NTHERM = 1 FOR PLOTS OF TEMPERATURE (DEG. C.) VS. WAVELENGTH 
C NTHERM = 0 FOR VOLTS VS. WAVELENGTH 
C NTHERM =-l FOR REFLECTANCE VS. WAVELENGTH 


MICHIGAN 
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C NTHERM --2 FOR RADIANCE VS. WAVELENGTH (RADIANCE UNITS = WATT/SQ.CM./ 
C CM./STERI X 10** ICST I 
C ICST = SCALING CONSTANT FOR THE PLOTS 

C LINE NO. = NS A > NS6, KS, NA, NR, KP, ID1, ID2 1615, 2A8) 

C LINE NO. .LE. 0 TO GET TO 'NO. OF FILES TO SKIP = • QUESTION 
C NC. OF FILES TO SKIP = -0 TO CHANGE CHANNEL SUBSET 
C NSf) = NO. OF LAST LINE OF INTEREST (NSB=-1 TO GO TO EOF) 

C NSK = NC. OF SCANLINES IN EACH STATISTICAL SET 

C NSK = 0 FOR STANCARO * I M PROVE * - TYPE RUN 

C ' I D 1 ' CAN BE STANOARD Afl IDENTIFICATION OR CAN BE THE TIME OF DAY 

C (FORMAT = XXXX.X) CORRESPONDING TO FIRST AREA 

C LEAVE 101, I D2 BLANK TO SIGNAL IRREGULARLY SHAPED AREAS, PUTTINC IN 
C ID ONLY WITH THE LAST LINE + POINT DESIGNATIONS 

C ID2 CAN BE STANDARD A8 IDENTIFICATION OR CAN BE 'DELTA TIME' (F6.3) 

C DELTA TIME = TIME INTERVAL BETWEEN AREAS (NEEDS TO BE ENTERED ONLY 
C ONCE AFTER PROGRAM IS CALLED UNLESS IT IS DESIRED TO CHANGE IT) 


DIMENSION A ( 1 9 ) , FLLI2A), LABIS), LABXI3), LBL!5,3>, PCILEI21), 

1 SC AN I 2 00 ) , T ( 12 ) 

COMMON DAT A ( 6000 ) 

COMMON FMTI19), ICI13), IT<9), LABEL(30>3), LINEI108), MINEUOO), 

1 PI 39 1, PP(3), QI39), QPI3) 

COMMON 8113,13), STDEVI13I, SUM ( 1 3 ) , NGOOD (13), FNI13.13) 

COMMON 80T ( 13), C0RRI13), ISCNI21,13), DATUM! 13), FLHII13), 

1 FLL0I13), IC0CEI13), ICP, NC, R, S, TDPI13), VI200.13) 

COMMON RS15), BANG2, DANG2, CC, NF, NR, MR, L80, ID1, ID2, 

1 BANG, DANC 

COMMON L90, KEY, NPTS, NUNES, NS A , NSB, KS, NA, NB, KP, IS, 

1 TITLEI12), TAG1I12), TAG2I12), NSS, NCHAN, KR, CONV, IPOS, 

2 IPACK, INT 

EQUIVALENCE (LINE, SCAN), ( T, TITLE) 

INTEGER DATA, FMT, Q, QP, R, REPLY, S, T, TAG1 , TAG2 , TITLE, U, W 

REAL 1 SON 

LOC (TEST = 70) 

CATAt A=1Q.» 8., 6., A., 3., 2., 1., .8, .6, .A, .3, .2, .1, 

1 .08, .06, .OA, .03, .02, .01) 

DATA I FLL = .Al, .A3, .A55, .A7, .A85, .50, .52, .55, . 58, 

1.63, .68, . 7A , .85) 

DAT A ( L A3 = ttH*.01* ,8H*.l * ,8H»1.0* ,BH*10.* ,8H*100* ) 

OATA(LA0X= 8H 12x , 8H 19X , HH 27X ) 

CATAILBL = 8H*2* 7X,8H*3* 6X,8H*A* 7X,BH*6* 6X,8H*8* AX, 

1 8H«2* 12X,8H*3* 9X,8H»A* 12X,8H»6* 8X,BH*8* 7X, 

2 8H«2« 16X,8H*3» 12X,8H*A* 16X,8H*6* 12X,8H»8* 9X ) 

DATA I PC I LE I 1 ) = 0., 5., 10., 15. ,20. ,25., 30. ,35. ,A0.,A5. ,50. ,55., 

1 60. ,65. ,70. ,75. ,80. ,85. ,90. ,95. ,100.) 

DATAISECOND = 1.) 

PARTMAP 


SLJ< 1) 


t>UU SL J I * } .INITIALIZE FOR NEXT REGION 

LOA(NC) SUB I 1 ) AJPI602I 

WRITE! A, 105 ) 

602 DO 6 10 JK= 1 , 19 
610 FMT(JK) = 8H 

CALL ZERCIB, FNI 169 ) ) 
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IP ~ 0 
IRREG - 0 
NFIRST = 0 
5L J ( 600 ) 

C SET UP TO CONSTRUCT LABELS FOR COMPUTER PLOTS 

1 CO 15 I L= 1 » 3 

LABEL ( 1 , 1 L ) = LAB! 1) 

DO 8 J L = 3 , 7 

* LABEL ( JLt [ L ) = LDL(JL-2,IL> 

LABE L ( 8 , I L ) = L AB ( 2 ) 

CD 4 JL = 10 » 14 

4 LABEL (JLt I L ) = LBL t JL-9 # I L ) 

LA BEL ( 15, IL ) = L A B ( 3) 

DU 5 JL= 1 7 i 21 

5 LABEL ( JL t I L J * LHL ! JL- 16 , I L ) 

L ABEL I 22 , I L > = LAB! 4) 

CO 6 J L- 24, 28 

6 LABEHJL, IL) = LBL ( J L-23 , I L ) 

CD 7 JL-2,30,7 

7 LABELS JL, IL) = L ABX ! ID 

15 L ABtL ! 29 , IL) = LABI 5) 

SL J4 t 600 ) - INITIALIZE FOR FIRST REGION 

DELTA = 0* 

K = 1 
L = 0 
MSA = 0 
N3 - 3 
NS R = 0 

101 FORMAT (615) 

102 FORMAT l 5E 15 .8) 

103 FORMAT (2013) ■ 

104 FORMAT ( A1 ) 

105 FORMAT I 1H1 ) 

10b FORMAT ( L2A8 ) 

109 FORMAT(lHO) 

118 FORMAT 1 1 Hi , 12A8, / 1 H , F4.2, 3H TO F6.2, I2H MICRON BAND /) 

121 FORMAT!/ 12X 1H0 10X 2H1. 10X 2H2. 10X 2H3* 10X 2H4. 10X 2H5. 

i 10X 2H6 • 1CX 2H7 • 11X 1 1HME AN ( STCE V > /) 

123 FORMAT < 1H 14A8) 

142 FORMAT ! F 8 • 1 ) 

144 FORMATUHO 2X 14HCHANNEL NO * = 10X 1317) 

145 FORMAT ! 8H LINES 15, 5H THRU 15, 8H EVERY 13, 10H, POINTS 14, 

1 5H THRU 15, 8F EVERY 13 *(* F6.2 * DEGREES TO * F6.2 ♦ DEGREES* 


LDAtSECONC) AJP( 12A) 

SECOND = 0, 

KR I TE ( 9 , 100 ) 

100 FORMAT ( 3H*#** JS1 ON ID TURN OFF EDITING* / 

1 3H**** J S 2 ON TO CONTINUE SEGMENTS ACROSS EOF* / 

2 3H*« ** JS3 ON TO ALLOW TITLE CHANGES* /) 
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1 2 A WKITEI9.107) 35B, 52B 

107 FORMAT ( * PUNCHED OUTPUT ON * 2R1) 

PE AD 19,1041 M2 

W = 0 

IF DATA CONVERSION ERROR CO TO 12 
CECOCEIl, 101, M2) W 

12 SL J4 ( 700 ) . SELECT INPUT MEDIUM 

NFIRST = 1 
WRITEIS, 108 ) 

108 FORMAT t • NC NSK ICST ICP NTHERM (515) * /) 
READ I R , 101 ) NC, NS 3 , ICST, ICP, NTHERM 

NC = M I NO I NC, 131 
ITHERM = NTHERM ♦ 3 
NTHERM = 0 

LDA(ITHERM) SUHI4) AJPU14) 

NTHERM = 1 


U CALL PR0CESS(DATA,1,0,N3,R ) 
IF I KEY .6Q.6 ) 12, 13 


13 

ISBAP = 7776 
LO A I 1 PCS) 

AJPKL + 2) 


LAC(O) 

S T A { 1 S BAD ) 


BACPT = FLOAT! ISBADI/CONV 


LDMNCHAN ) 
NC = 1 
NC i * 1 
NC2 = 1 

l N A ( -2 ) AJPZI94) 

94 

LDA(NC) 

AJPK 11) 


NC = M l NO ( 1 3* 

NCHAN) 

11 

NC3 = MAX 0 ( 1, 

NC/2) 


NC4 = M I NO < NC • 

NC 3 + 1 ) 


NCI = IC0DEINC3) 
NC2 = ICOCE(NCA) 


C NCI ♦ NC2 ARE USED IN EDITING,.. IF THE SIGNALS IN BOTH OF THESE 
C USUALLY STRONG CHANNELS ARE 'RAD', POINT WILL BE DISCARDED 

MR IT E t S , 1 10 ) NC 

110 FORMAT!* SPECTROMETER CHANNELS (* I 2* I 3 1 = *) 

READ I R , 103 ) I 1C I I ) , I = 1, NC) 

FLLO(l) = 0. 

DO 10 1=1, NC 
ICGDEtU = I 
LDAKICI AJP(IO) 

1TMP = ( 1 6* I T MP ) / 1 5 
IC I I ) * 118*1 TMP ) / 1 7 
FLLOl I ) = FLLUCIII) 

FLHII 1) = FLLIICm+1 ) 

10 CONTINUE 

LOA(NC) SU8INCHAN) AJPtlTA) 

MRITE(S,lll) NC 

111 FORMAT ( » SUBSET OF CHANN EL S ( * I 2* I 3 ) = *) 

RE AD (R, 103) < I CODE t I ), 1 = 1, NC) 

13A SLJ4I700). SELECT INPUT MEDIUM 


36 




2 CALL PROCESS (DATA, 2) 

I N T = 1 
ANGA = BANG 2 

ANGB * ANGA + ( FL0AT1NPTS 1-1. )»OANG2 
LOA(KEY) SUB ( 5 ) A J P ( 12) 

IF(ID1 .EC. 1H ) 21, 20 

20 IF DATA CONVERSION ERROR GO TO 21 
CECOCEI 1,142, 101) TEMP 

CECOCEI 1,142, 102) DEL 
LDA(CEL) A JP ( 2 1 ) 

IU2 = 1H 
DELTA = CEL 

21 NS K = NS3 
AJP1I22I 
NS F = NSB 

NSK = NSB - NS A 1 

22 IF! I Cl .EQ. 1H ) 23, 24 

23 IRREC = 1 

WR 1 TE ( 4 , 145 ) NS A , NSF , KS, NA, N9» KP, ANGA, ANGB 
SLJ125) 


24 

IRRFG = 0 


25 

SLJ2I26) 

SUC27) 

26 

MSA = NSR 
NS R - 0 


27 

NSF - NSA + 
MSA = 0 
EDIT = 10, 
suit 29) 

EDIT = 0. 

NSK - 1 - MSA 

29 

NFIRST = 1 



SL J3< 39) 

SL J 1 37 ) 

39 

WRITEI9* 117) 


117 

FORMAT { * NEW 

TITLE * /) 

40 

READ ( 9 f 106) 

IT 


iFiirm *eg - 

3HYES ) 40, 41 

41 

1F( I T C 1 > .ECU 

2 FNO ) 37, 42 

42 

CO 43 1 = 1,9 


43 

T [ TLE t I ) - 

IT (I) 


37 IFINC .EC. 1 .AND. NFIRST .GT, 0) 44, 45 

44 RRITF(4, 118 > TITLE, FLLO(l), FLHK1) 

SCALE = 12. 

F THERM = 1. 

OFFSET = 1. 

45 MGCOD = < ( NSF-NSAI/KS + l )• ( INB-NA )/KP4l ) 


3 CALL PROCESS (DATA, 3) 

LL = NCHAN* I NA- 1 ) 

CO 49 1SS=NA,NH,KP 

LOA(IP) SUB ( 200 ) A JP 3(47) 
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CO 46 J=1,NC 

L = LL + l COO E C J ) 

LDA4!DATA> 4JP1I46A) 

ID A | l SB AD I STA4IDATAI 

46 A DATUM! J ) = FLOAT (DAT AIL ! I/CONV 
46 CONTINUE 


LL = LL*NCHAN»KP 

SL J4 t 200 ) . TEST POINTS AND CUMULATE TOTALS 

SL J ( 49 I 


47 IP = IP ♦ l 

CO 48 J-liNC 

L = LL ♦ I CODE ( J ) 

LDA4 ( DAT A I AJP1I48A) 

LO A ( I S BAD I ST A4 I DATA ) 

48 A V < I P , J ) = FLGATIOATAIL D/CONV 

48 CONTINUE 

LL = LL+NCFi AN*KP 
IFIIP .EC. 2001 51, 49 

51 SLJ4I300). GET BOUNDS FOR COLLECTED POINTS 

49 CONTINUE 

IF(NSF-IS) 52, 52, 50 

50 I F ( KEY-3 ) 3, 52, 90 

52 NGCOC = NGCOD + MGOOO 

LOA(IRREG) AJP 1 ( 2 ) 

LDA(IP) SUB 12001 

SLJ4I300). 

53 DO 54 1*1, NC 

54 SUM! I ) = SUM! tl/FN! 1,1) 

CO 55 1=1, NC 


CO 55 J= I , NC 



B ( 1 » J ) « B( 

If J > /FN ( I , J) 

- SUM! I ) 

•SUM! J ) 

55 

ft! J, I ) = H! 

It J) 



LOAINC) 

SUB ( 1) 

AJP1I 

60) 


LOA(NFIRST) 

A J P l 56 ) 




LOA { NTHERM ) 

A J P < 56 ) 




HRITE(4, 120 ) 

120 FORMAT!/ 10X 3H-4Q 9X 3H-30 9X 3H-20 9X 3H-10 10X 1H0 10X 2M1 n 
1 10X 2H20 10X 2H30 10X 3H40 1 1HMEAN ( STDE V ) /) 

OFFSET = 61. 

FTHERM * 10. 

56 P ( 2 ) = SUM { 1 ) 

B ( 1,11 = AM AX 1 ( B I 1 , 1 } , 0.) 

STCEV(l) * SORT ( B1 1, 1 ) ) 

P!l) = PI 2 ) - STOEV(l) 

P! 3) = P( 2) * STOEVI 1 ) 

FSUM = SUM! 1 ) »F THE RM 
FSTDV = STCEVtl )*FTHERM 


AJP2! 53) 

GET SOUNDS FOR COLLECTED POINTS 
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CO 57 J= 1 , 3 

Q I J ) 3 P( J ) *SCALE + OFFSET 
blJI = MINOIQI J) , 108) 

57 Ql J) = MAXO(CHJ) , 1) 

CO 58 J=1 , 108 

58 LINEIJ) = 1H 

UNEIOIllI = 1 H ( 

LINE I 0(3) ) = 1H) 

L I NE ( Q I 2 I I = 1HX 

LOAINfHERM) SUB I 1 ) AJPI59) 

LDA(NFIRST) AJP ( 59 ) 

WRITCI4.121 ) 

59 HR I TE ( 4 , 122 ) 101, ILINEIJ), J=l,96)» FSUM, FSTDV 

122 FORMAT ( ItHO IX A8, 2X 96A1, F6.2, *(• F4.2 •)* /) 

SLJI65) 

60 LDA(ICP) A J P ( 6 1 ) 

LOA(NFIRST) AJP(61) 

HRITE<5,118) T, FLLOIICP). FLHItlCP) 

N5 = 5 

SLJ4UOO). IDENTIFY UNITS OF COMPUTER PLOTS 

61 WRITEI4.123) T, 101, ID2 

WRITEI4, 145) NSA, NSF, KS, NA, N8 , KP, ANGA, ANGB 
WRITEC4, 144) (ICODE(I), 1=1, NC) 

WRI TE I 4 , 124 ) INGOODC I), 1 = 1, NC) 

124 FORMAT! 12HOGOOC POINTS 17X 1317 /) 

HR I TE I 4 , 125 | I SUM I 1 ) , 1 = 1, NC) 

125 FORMAT ( 12H0MEAN VECTOR 17X 13F7.3) 

SINC = 0. 

CO 62 1=1, NC 

!F(B(!,I) ,LT. 0.) 63, 64 

63 SING = 1. 

BU.II = 0. 

STCEVI 11=0. 

SL J ( 62 ) 

64 STCEVI I ) = SORT I B( I, I) ) 

62 CUNTINUE 

HR I TE ( 4 , 126 ) (STDEVIt), 1 = 1, NC) 

126 FORMAT ( 20H0 ST AN CARD DEVIATION 9X 13F7.3) 

WRITEI4, 127) 

127 FORMAT 1 19H0CCRRELAT I ON MATRIX ) 

CO 99 1=1, NC 

CALL ZERO (COR R( 1), CORRt I ) ) 

CORRI 11= 1. 

LDAl(STOEV) A J P ( 67 ) 

II = I-i 

CO 66 J = 1 , I 1 

LDA2ISTDEVI AJPI67) 

66 CORRIJ) = BIJ,I )/STOEV( I l/STDEVIJ) 
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67 LDA3 I FLLC1 ) AJPI97) 

WK I TE I A , 1 28 ) FLLO(I), FLHIII), ICORRIJ), J = l,l) 

128 FORMAT ( 1 HO IX FA.2»AH TO F5.2, 1 AH MICRON BAND 13F7.3) 

SLJI99) 

97 WRITEIA, 1A6 I ICCDEII), ICDRR(J), J=1,I) 

1A6 FORMAT I 1 HO 2AX 12, 2X 13F7.3) 

99 CONTINUE 

LDAINGCOD) SUBINCI AJP2I68). TEST FOR SINGULAR MATRIX 

WRITEIA, 129) NGOOD(O), NC 

129 F0RMATI1H0 12, * POINTS AND* 13, • CHANNELS, SO SINGULAR * 

1 *CCV AR I ANCE MATRIX *) 

SL J< 69 ) 

68 LDA(SING) AJ P ( 70 ) 

WRITEIA, 130) 

130 FORMAT!/ » SINGULAR MATRIX *) 

SLJI69) 

70 CALL POSCEF I B I 0, 1 ) » NC, 13, 1, 1H0, SUMIOJ) 

69 WRITE I A , 131 I 

131 FORMAT ( 1H 1 *MEAN VECTOR AND COVARIANCE MATRIX * /) 

WRITEIA, 102 ) ( SUM ( I } , 1 = 1, NC), 1IB(I,J), J=I,NC), 1 = 1, NC) 

WRITE I A, 109) 

WRITEIA, 1AA ) IICODEIJ), J - 1, NC ) 

WRITEIA, 113) IBOTI I ), 1=1, NC) 

113 FORMAT I 1 HO 8X»LCWER BOUND EDIT*2X 13F7.3/) 

WRITEIA, 112) 

112 FORMAT I 1H0 11X*PERCENTILE*/12X» */) 

00 1 1 A J J= 1 ,21 

WRITEIA, 1 1 A ) PCILEIJJ), II SCN 1 J J , J ) , J = 1 , NC > 

11A FORMAT [ 11X, F8.1, 8X 13F7.3) 

WRITEIA, 115) I TOP I I ) , 1 = 1, NC) 

115 FORMATUHO aXtUPPER BOUNO E 0 1 T* 2X 13 F 7 . 3 /I 
WRITEIA, 105) 

WRITEIA, 123) TITLE, 101, ID2 
WRITEIA, 1321 

132 FORMAT! 1H0 33X*SEMIL0G PLOT OF MEAN +-1 STANDARD DEVIATION*) 

LDA3IFLL0) AJPI71) 

WRITEIA, 133 ) 

FQRMAU1H + 10AX 15HW AVELENGTH BAND / 108X 9HIMICR0NS) ) 

71 WRITEIA, 13A) 

13A FORMAT! 1H0 108X 7HCHANNEL ) 

72 P2 = 0. 

PI = 100. 

CO 73 1=1, NC 

P2 = AMAX 1 1 SUM ( I ) , P2 ) 

73 PI = AM IN 1 1 SUM! I ) , PI ) 

P2 = 1.1*P2 
PI * 0 .9* P 1 
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P2 = AKIN H 100. < P 2 ) 

PI = AKAXU0.01, PI) 

FSCALE = ALOGIOI P2/P1 I 
NSCALE = A. /FSCALE 
NSCALE = M I NO I A » NSCALE) 

NSCALE = MAXOI2, NSCALE) 

SCALE = 2 5 »NSC AL E 
FMID = SORT ( P2 »P 1 ) 

FSCALE = 50. /SCALE 
PA = FMID/UO.R+FSCALE) 

PA = AKAXKPA, 0.02) 

EN1AI19I LDA(PA) 

+THSAIA+1) SL J ( 92 ) 

JAC = L + l 

I AC = ! ( 18-JAC)*7-1 ) / 6 

I AD = HAXOI I AC, -1) 

SHIFT » I, - AL0G10IAI JAD+l ) )*SCALE 

00 60 1*1, NC 

I I = I COD El I ) 

P I 2 I * SUMII) 

PI 1) * PI 2) - STDEVI I ) 

P t 3 ) = P ( 2 ) ♦ STDEVI I ) 

LDA(NFIRST) A J P ( 78 ) 

P < A ) * 0.01 

CO 7 A M*5» 1 3 

7A PlM) = PI M-l) ♦ 0.01 

CO 75 M= 1 A, 22 

75 PIP) * PI M-l ) + 0. 1 

CO 76 M*23, 31 

76 PtM) = PIM-l) + 1. 

CO 77 M= 32 , 36 

77 P(P) = PI M-l) +10. 

78 CO 81 J= 1 , 36 

L0A2IP) AJPI79) AJP3I79) 

QIJ) * ALOGlOlPtJ) >*SCALE ♦ SHIFT 
Cl J) = Ml NO I Q( J ) , 108) 

Cl J) = MAXOIQI J ) , 1 ) 

SLJI81) 

79 GUI * 1 

81 CONTINUE 

DO 82 J= 1 * 1 08 

82 L1NEIJ) * 1H 

CO 83 M*A, 36 

83 LINE(QIM) ) = 1H. 

LINEI108) * 1H 
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CO 36 I PG=Q( 1) , c 1 3) 

36 LINEUPQ) = 1H- 

LINE!Q(1) ) = 1HI 
l I NE ( QI 3 ) > = 1HI 
L 1 NE ! Q 1 2 ) ) = 1H» 

L0A3IFLL0) AJPIB4) 

WR I T C ! 4 * 1 35 ) LINE, FLLO(I), FLHKII 

135 FORMAT I / 1H0 108A1, F4.2, 2H - F5.2/J 

LOA(ICP) AJPOO) SLJIB5) 

84 WR IT E 1 4 , 1 36 I LINE, II 

136 FORMAK/1HO 10HA1, 2X 12 /} 

85 10 A ( I I ) SimUCP) AJP1I80) 

PP12) = P ( 2 ) 

PP(l) = P(l) 

PP(3) = P( 3) 

CO 86 M* 1 , 3 

CP ( M J = PP!M)*12. + 1. 

GP(M) = M I NO < CP!M), 100) 

86 CP ( M ) = MAXOI CP 1 M) , 1 ) 

CO 87 J=l,100 

87 MINE(J) = 1H 

LDA(NFIRST) AJPI88) 

WR I TE t 5 , 121 ) 

88 CO 38 IPOQPtl) ,CP( 3 ) 

38 MINE! IPO = 1H- 

MINE(CPIl)) = 1HI 
MINE ( QP( 3 ) ) = 1HI 
M I Nfc < QP ( 2 ) ) = 1H* 

KR ITE 1 5 * 122 ) ID1, (MINEIJ), J=l,96), SUM(I), STDEVtl) 

80 CONTINUE 


MJM = 28/ NSC AL E + 2 
FMT(l) = 8H< IX 

CO 89 JM=2, MJM 

89 FMT(JM) = LABEL! JM+IAO, NSCALE-1) 

FMTIMJM+1) = 8H ) 

WRITE!4,FMT ) 

N5 = 4 

SLJ4I800). IDENTIFY UNITS OF COMPUTER PLOTS 

LDA { W ) AJP ( 65 I 

WR I TE ( W , 140 ) ITITLE(I), 1 = 1,9), ID1, 578 
140 FORMAT! 8A8, A6, 2X, A8, IX, Rl) 

ENCOCE(10, 102, DATA) I SUM! I),I = l,NC),((B!I,J),J = I,NC),I = l t NC) 
NCARD = INC + NC*(NC+l)/2 + 4)/5 
J = 0 
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CO 155 1 = 1, NC AR C 




WR IT E ( W , 141 1 (DATA! L ) , L = J, 

J+9), I, 57B 

141 

FORMAT ( 9A8, A3, 3X, 12, IX, 

R 1 ) 

155 

J = J ♦ 10 



65 

LDA { NSF ) 

SUB [ NSB ) 

AJP1I156) 

161 

SL J4 ( 600 ) • 


RE-INITIALIZE FOR NEXT REGION 


SL J ( 2 ) , 


CALL PROCESS! DATA, 2) 

156 

NSA = NSF ♦ KS 




LDA ( NSB ) 

SUB I NSF) 

SUB ( NSK ) AJP2I157) 


NSR = NSB - NSF 




NSF = NSB 




SLJ< 158 ) 




157 NSF = NSF + NSK 

158 LDA ( DELTA ) AJPI160) 

CECOCEtl, 142,101) X I D 
XIO = X ID + DELTA 
159 ENC0CE(1,142,ID1) XIO 

CECOCEll, 143,101) IH 
143 FORMAT | 4X , II, 3X ) 

LDAUH) SUB 1 6 ) AJPK160) 

X I D * X I D ♦ 40. 

SLJ( 159) 

160 SL J4 (600 ) . RE-INITIALIZE FOR NEXT REGION 

IF(NSF-NSA) 2,2,45 


90 NS B = IS - KS 
NSF * NSB 

NSR = NSF - NSA ♦ l 

MGCOD = < ( NS6-NSAJ/KS+1 )• ( (N8-NA J/KP+l) 
LOA(NSR) SUBI2I 

A JP3 ( 1611 SLJI52) 


200 SL J I • ) . TEST EACH POINT. IF IT PASSES THE TESTS,' CUMULATE SUMS 
I F ( DATUM (NCI) . EQ ■ BADPT .AND. DATUM ( NC2 ) .EQ. BADPT) 200, 201 

201 SLJK203) 

DO 210 J = 1 * NC 

I F ( DATUM ( J ) . GT . TOP(J) .OR. DATUM ( J ) .IT. BDT(J)) 200, 210 

210 CONTINUE 

203 DO 230 1=1, NC 

LDA1 ( DATUM ) SUBIBADPT) AJPI230) 

SUM ( I ) = SUM ( I ) + DATUM I 1 ) 

NGCOD { I ) = NGCOOI I ) + 1 

DO 220 J= I , NC 

L0A2 ( DATUM ) SUBIBADPT) AJPI220) 

8 1 I , J ) = 0 1 I , J ) + DATUM! I >#DATUM( J) 

FN( I , J) = FNI1.J) + 1. 

220 CONTINUE 

230 CONTINUE 
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SL J ( 200 ) 

300 SL J ( • ) • GET BOUNDS FOR COLLECTED POINTS, THEN CUMULATE SUMS 
SLJH311) 

CO 310 J=l,NC 

CALL MEDI ANIVIO, J ), SCAN, IP, CENTER, QDEV ) 

ISCNll.J) = SCAM 1 ) 

DO 330 J J = 1 , 20 

NUMB = J J • ! IP-1 1/20 + 1 
330 I SC N ! J J + 1 , J ) = SCAN ( NUM B ) 

ODEV - AMAX1IC0EV, 0.867/C0NV, EDIT) 

TOP ( J ) = CENTER * 5.77*QDEV 

310 BOTU) = CENTER - 5.77‘QCEV 

C .001 LEGITIMATE PTS. EXCLUDED. VARIANCE DECREASE .LT. 1 0/0 

311 CO 320 N= 1 , IP 

CO 315 J* 1 , NC 
315 DATUM ( J ) = VIK.J) 

SLJ4I200). TEST POINTS AND CUMULATE SUMS 

320 CONTINUE 

SL J ( 300 ) 

700 SL J ( • ) . SELECT INPUT MECtUM 

WRITE19, 7011 

701 FORMAT!* INPUT MEDIUM* •) 

RE AO ( 9 , 104 ) REPLY 

LOA(REPLY) AL S ! 6 ) 

STAIR). R= INPUT MEDIUM 

S = 9 

I F ( REPLY .NE. 1HT ) 702, 700 

702 LIU1ITEST) IJP1I700) .IF BATCH, SAVE TYPEWRITER OUTPUT 

S = 23 

SL J ( 700 ) 

BOO SL J ( * I . IDENTIFY UNITS OF COMPUTER PLOT 

GO TO 1153, 152, 151), ITHERM 

151 WRITEIN5.137) 

137 FORMAT!// 1H9 56X *VOLTS* ) 

SL J I BOO ) 

152 WRITEIN5.138) ICST 

138 FORMAT!/ 1H9 67X, 12 / 1H9 52X ‘REFLECTANCE X 10 * I 

SL J ( BOO ) 

153 WRITEIN5.139I ICST 

139 FORMAT!/ 1H9 77X, 12 / 1H9 42X ‘RADIANCE ! WA TT/ SQ . CM. /CM. / STER ) X 
110 * ) 

SL Jl 800 J 

92 CALL STOP! 7HBAD RUN I 
END t PAGE 
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FUNCTION PCSDEFIA, NN» IDA, ISW, CC, MEAN, N TAPE I 
C USES 2*N»«2 ♦ N ♦ 3 SPACES IN ERASABLE 
DIMENSION A I I DA , 1 DA ) , MEAN! IDA), Ell) 

REAL MEAN 

INTEGER CC, L, N, 0, Q, 0 
CA T A ( L = A ) 

LOC I ERAS = 63) 

LINl I I , JJ) = N»t JJ-1 I t 1 1 

N * NN 

IF tN.ee. 1 } 25, 2 
2 N 1 = N +■ 1 

N2 * N*N ♦ l 
ENA I E I SUB ( ERAS ) 

SCM(-OB) ST A I B ) 

ADD I N 1 ) STAIC) 

ADDIN2) STAID) 

ENA(NTAPE) AJP(l) 

L = NT APE 

1 CALL MXMCVIO, A, BID), N, N, IDA, N) 

CALL JACCBIN(EIC) ,N,N,E(B), EIO), -L) 

WRITER, 105) CC 

105 FORMAT I Al, *E I GEN VALUES IYVAR) YMU COEFFICIENTS OF THE 

1X*S * /) 

POSDEF = 0. 

DO 24 I = 1,N 

VAR = EID + L IN I I , I ) ) 

IFIVAR) 5, 5, 6 

5 POSDEF * 1. 

6 YMU = 0. 

IF(ISW) 10, 15, 10 
10 CO 12 J = 1 ,N 

12 YMU * YMU + E(Q + L IN ( J, I ) ) *MEAN I J ) 

15 WRITE(L,115) VAR, YMU.IEIC ♦ LINIJ,IM, J = 1,N) 

115 FORMAT I 1H 2E13.5, 2X 13F7.4 / < 29X 13F7.AJ) 

2A CONTINUE 
25 RETURN 

END 


MACHINE A2F 3 ICHAN, SMOOTH, A, B , C , D , E, F , G , H ) 

LOC ( Z=0,P0 INT=30,ERP=24,MTR=25,RNX=11, RDB=100,WRB=300,CNT=60) 
RSV ( A 1 = 0 ) 


SIMPLE MESSAGES 


CON (MES1=2062612363702040B*N1 
CON (MES 2=2 06 10 266206 15 167B, N2 
CON (MES4=2023712343653732B,N3 
CON (MES6*2020204461453020B| 
CON I MES COR E = 2223 A6 5 161676520B 
CON IW0NKITY=20264645427123306 
CON IMA =A72A23202361A765B 

CON ( M6 =6A30730000000000B 

CON {DONE =6AA6A56520267123B 

CON (WARN =20A6632361A32061B 


=20A5A62061216A0OB)» 

=20655151A65100008)» 

= 00000000000000006) . 

,M2=A371AA7123000000B 

,M3=1A52733220A62A236 

,M5=20A5A62320516561B 

) 

, M 7=70 20 2361 A7653700B 
,M8=7 16722 73000000008 


TAPE SELECT CODES 


A2F30000 
A2F30010 
A2F30020 
A2F30030 
A2F3O0A0 
A2F30050 
A2F30060 
A2F30070 
BATCH - NO A/D A2F30080 
A2F ARG ERROR A2F30090 
TITLE A2F30100 

MANY A2F301 10 

). STORAGE LI MI TA2F30120 
(.OUTPUT TAPE A2F30130 
) . NOT READY A2F301A0 
A2F30150 
(.DONE WITH TAPA2F30 160 
). OCTAL ARGS. A2F30170 
A2F30180 
A2F30190 
A2F30200 
A2F30210 
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CON ( T0=0 ) 

CON (T 1=014000000004201 IB* T2=024G0CC00Q042O21B J 
CON (T 3 = 034 0000 00004 203 IB* T 4 = 0440000000042 0418 ) 
CUN (T5*054U00000004 30 11B,T 6 = 064000000004 302 IB) 
CON (T 7 = 074 000000 004 30 3 IB, 18=1040 00000004 304 IB ) 


-DUMMY 
.1 2 
• 3 4 
. 5 6 
.7 8 


L0C 

CON 

CON 

CON 

CON 

CON 

CON 

CON 

CON 

CON 

CON 

CON 

CON 

CON 

CON 


( NUMST I X= 14 ) 

ICRRGRS=0, SHORT = 0* BURNED=0, 


LO ST-0 t LONG =0 * STRETCHED) 


( SC ANS= 0 > L I NES = 0 * SIZE=0, NS$=0, NCHAN=0, SM0QTH=0> 


A2F30220 
A2F302 30 
A2F30Z40 
A2F30250 
A2F3026O 
A2F30270 
A2F30280 
A2F30290 


* MASKS 

* 

f * 







A2F30300 
A2F30310 
A2F303 20 

CON 

( MAS K 1 1 =0 

777 

000 

777 

000 

777 

B ) 

A2F30330 

CON 

( MASK 12 =0 

377 

000 

377 

000 

377 

B ) 

A 2F 30340 

CON 

( MASK21 =0 

000 

777 

000 

777 

000 

6) 

A2F30350 

CON 

( MASK22 =0 

000 

377 

000 

377 

000 

B ) 

A 2F 30360 

CON 

( MASK3= 4 

003 

777 

777 

777 

777 

B ) 

A2F3037Q 

CON 

( M AS K0= 4 

003 

774 

003 

774 

003 

B ) 

A2F30380 

CON 

<MASKE= 7 

774 

003 

774 

003 

777 

B ) 

A 2F 30390 

CON 

(MASK =0 

400 

400 

400 

400 

400 

6) 

A2F30400 

A2F30410 

CON 

(MASKP1 - 

000 

000 

000 

000 

777 

OB ) 

A2F30420 

CON 

( MASKP2= 

000 

000 

000 

77? 

OOC 

Ob) 

A2F30430 

CON 

l MASKP3 = 

000 

000 

777 

000 

000 

OB) 

A2F30440 

CON 

(MASKP4= 

oco 

777 

000 

000 

000 

OB) 

A2F30450 

CON 

( M A S K P 5 = 

777 

000 

000 

000 

000 

06 ) 

A2F30460 

A2F30470 

CON 

<MASK1=0, 

M A S K2= 0 ) 





A2F30480 









A2F30490 




»»»*« 

■* 




A2F30500 

* JCB STATISTICS MESSAGES 

• 




A2F30510 

»«*#»*+««« 

!■#**###*##< 

nt # 

m- * * « 

* 




A2F30520 


( JJ=0, 
l J3-0 $ 
( J6=0, 
< J9 = 0* 
l JB=0» 
< JE = 0, 
( JH = 0, 
l JL=Q , 
{ J0=0, 
t JR=0, 


J1=236147652G6551516, 
J4=2C2 02 34646202270B, 
J7=2020612320222361B* 
J0=2 04346222 3000 000B, 
JC=2C23517144446564B, 
JF =2022235 165236 370B, 


A2F30530 
A2F30540 
A2F3O550 
A2F30560 
A2F3O570 

J 2=46 5 122 0000 000 OOQB) . TAPE ERRORSA2F30580 


TOO SHQRTA2F30590 
AT START A2F30600 


J5=4651230000000000B) 
J8=5123000000000000B) 

JA=C0000000000000006 ) - LOST A2F30610 

JD=U0000000000000006) . TRIMMED A2F30620 
J G= 6 564-OOOOOOOOOOOOB ) • STRETCHED A2F3063Q 
J I=22636l4522000000B, J K = 00 00 00000 OOOOOOOB ) . SCANS A2F 30640 

JM=4 37 1456522000000B, J N = 0000 00000 0000 00 OB ) .LINES A2F30650 

JP=26465l6422214371B, JQ=4565 0000000000006 > .WORDS/LINE A2F30660 
JS=2261444743652221B, J T = 4 37145 65 0 00000006 ) . SAMPLES/L I NA2F3067O 
t JU = 0, JV=6 3 706 1454 56 54 32 OB , JW =6 46 1236100 00 000 OB ) .CHANNEL DA TA2F 30680 
(JX=0, JY=2 346200 1202244466* JZ=4623 70000 00000 OOB ) . TD 1 SM00THA2F3069O 


A2F30700 

******** ,*,*»,»* A2F30710 

* UlHtR JUNK * A2F30720 

A2F30730 

A2F30740 

CON( CNDREEL =17654051656543175 ) A2F30750 

CdN l SPACES =202020202020202 OB) A2F30760 

CON ( FIRST =0000000 1 777 77777B) A2F3077Q 

HOUSMOTHD = SMOC THED ) A2F30780 

A2F30790 

***** A2F 30800 

* BEGIN MAIN PROGRAM * A2F30810 

************************ A2F30820 

A2F30830 

LENT SLJ (*> SLJ (l+ll) .EXIT/ENTRY A2F30840 

OARG BSS { 1 ) .CHANNELS A2F30650 
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9 ARG BSS U) 

1 ARG BSS (8) 

L l LI (POINT) 
IJP1U+1) 
LOGKZ + 2) 

ENI 1 ( ERP ) 

BSS (1) 

ENA 18) 

THS I ( 1 ARG ) 
1AER ENA (0) 

STA1 (1ARG) 
ENA IMES2) 
ENA (0) 

THS1 ( 1 ARG ) 
SIU1 ( 2E0T) 

2 ARG LDAl(lARG) 

AJP ( 1 AER ) 
SILKL+U 
STC1 ( 9S EL ) 
EQS311ARG) 

SLJ ( 1 AER ) 

ENA (0) 

I JPK2ARG) 

ST A1 (ERRORS ) 
LDA (OARG) 

THS (16) 
AJP3UAER) 

SAL ( 4 1 A } 

SAU (6IA) 

SAU (2TST) 

SAU ( 2 1 A ) 

SAL I1IA) 

ENA (A 1 + 1) 

SCM I-OB) 

LDA { 9ARG ) 

THS (16) 

INA ( “3 ) 

INA (2) 

EN I It 2 ) 
THSUOARG) 

ENA (WARN) 
ENIKO) 

ENA (0) 

SCA6 ( 2 ) 

I JP6IL+1) 

I JP6IL + LJ 
SIU1 ( 100S) 
LDQ1 (MASK 11 > 
LUGl ( MASK21 ) 

S I L6 { 200 ) 

STA (IKR) 

ENQ6 ( — 1 74 16 B ) 
CJP61ERP) 

5AL (Z + 5) 

SCM (-06) 

LDA (Z + 7) 

SAU UOUT) 

ENA (LENT) 
MEQ7 ( CNT *7 ) 
INIl(l) 

ENA (Al) 

SIL ( LFLG ) 

SIL ( 1 EOT ) 

SAL (7NXT) 


ENA l ME SI ) 

SLJ IL+4) 

STQ (L + 2) 
SILUL+1) 

ENI1I8) 

SLJ ( L + 4 ) 

ENI 1 ( 7 ) 

I JP1IL) 

SLJ4( ERP) 

ENI 1 ( B ) 

SLJ (1AER) 

ENI2< TO) 

SAL (L+2) 

A JP 3 ( 1 A ER ) 

LDQ2 ( * ) 

EN I 3 ( * ) 

SLJ (L+2) 

STA1( I ARG ) 

ENI 1 ( NUMST t X- 1 ) 
IJPl(L) 

AJP (LAER) 

SLJ ( 1AER ) 

STA (NCHAN) 

ALS m 
INA (1) 

SCM ( —0 B ) 

LAC INCHAN) 

SAU ( 3 I A ) 

ADD (NCHAN) 

SAU 11AGN) 

STA (SMOOTH) 

SLJ ( 1 A ER ) 

A JP 3 ( 1 A ER ) 

SAU ( 10 IE ) 

ENA (7) 

SLJ (L+2) 

SLJ4{ MTR ) 

ENQ (L) 

DVI (SMOOTH) 

$CL ( MASK3 ) 

ENI 1 ( 1 ) 

EN I 6 ( 1 ) 

SIU1( LOIS) 

STQ (MASK1) 

STQ IMASK2) 

S I L 6 ( 20 OM ) 

EN16( Al ) 

ENA (MESCORE) 
ENA6( 1200) 

SAL ( OSTX ) 

SAL ( 2K ) 

SAU (9E0T) 

ENQ ( 37777B ) 
LILllCNT+5) 

SLJ (L+3) 
SlLKCNT+5) 

SAL (CNT+3) 

SIU ( 8NXT ) 

ENA ( 37 777B ) 

SAU ( 9NXT ) 


SMOOTHING FACTOR 
TAPES 

TEST FOR NEW RESIDENT 
GO AHEAD IF OLD RESIDENT 


EXIT IF IN BATCH 

OK. DO ARGUMENTS 

ERROR IF ARGUMENT LARGER THAN 8 

ARGUMENT ERROR 

CLEAR FOR NEXT CALL 

TYPE MESSAGE/EXIT 

GET NUMBER OF ARGUMENTS 

SAVE COUNT, GET ADDRESS OF TABLE 
GET ARGUMENT 

ERROR IF ZERO OR NEGATIVE 
SET SELECT CODE 
SAVE COUNT 

SEARCH FOR DUPLICATE ARGUMENTS 
ERROR IF DUPLICATE ARG FOUND 
CLEAR ARGUMENT 
LOOP FOR MORE TAPES 
CLEAR JOB STATISTICS 
ERROR IF NCHAN , EQ. 0 
ERROR IF NCHAN .GT. 16 
IF NCHAN .GT.-O, STORE IT 


STORE SMOOTH FACTOR 
ERROR IF SMOOTH .GT. 16 
ERROR IF SMOOTH .LE. 2 
SET UP SMOOTH LOOP LIMIT 

TEST FOR SMOOTH OR CHAN OVER 7 
ISSUE WARNING MESSAGE 


CRASH IF TOO NEAR END OF CORE 
STORE END OF INITIAL READ BUFF 


SAVE INTERRUPT ROUTINE ADDRESS 

FIND A2F IN TABLE 
REMOVE REMAINING ENTRIES 
SET END OF A 2F 
STORE ZERO 

STORE 37777 


A2F3O06O 

A2F30870 

A2F30880 

A2F3089Q 

A2F30900 

A2F309 10 

A2F30920 

A2F30930 

A2F30940 

A2F30950 

A2F30960 

A2F30970 

A 2F30980 

A2F30990 

A2F31000 

A2F31010 

A2F31020 

A2F31030 

A2F31040 

A2F31050 

A2F31060 

A2F31070 

A2F31080 

A2F31090 

A 2 F 3 1 1 00 

A2F31U0 

A2F3U20 

A2F31 1 30 

A2F31140 

A2F31150 

A2F31160 

A2F31170 

A2F31180 

A2F31190 

A2F31200 

A2F31210 

A2F31220 

A2F31230 

A2F31240 

A2F31250 

A2F31260 

A2F31270 

A2F312B0 

A2F31290 

A2F31300 

A2F31310 

A2F31320 

A2F31330 

A2F31340 

A2F31350 

A2F31360 

A2F31370 

A2F313B0 

A2F31390 

A2F31400 

A2F31410 

A2F31420 

A2F31430 

A2F31440 

A 2F 31 450 

A2F31460 

A2F31470 

A2F31480 

A2F31490 
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SAL (2FLG) 

SAL ( OF LG ) 

•STORE NON-ZERO 

A2F31500 

ENA (WRft) 

SAU (1WRT) 

•MORE ADDRESSES 

A2F3151Q 

ENA (WR6+54) 

SAU (OWRT) 


A2F31520 

SAL IZ + 4) 

ENUI Al > 

* 

A2F31530 

ENQ (0) 

STQKZ) 

.CLEAR CORE 

A2F31540 

I SK 1 ( 37777B ) 

SLJ (L-l) 

9 

A2F31550 

ENA ( 2WR T ) 

SAU (OTAPJ 

• 

A2F3156D 

ENA (OENC) 

SAL ( IT YP ) 

« 

A2F31570 

ENA ( 0 1 N ) 

SAU (0A2C) 

* 

A2F31580 

EXF uoitn 

EXF ( 2000B ) 

•CLEAR FAULT INT., STOP CLOCK 

A2F31590 

A2F31600 


a • 


A 2F 3 16 10 

* FIRST TAPE RECORD 

* 


A2F31620 


« * 


A2F31630 

A2F31640 

ENA (MESA) 

SL J 4 ( MTR ) 

• TYPE * T I TLE* 

A2F31650 

ENA (9) 

SL J 4 ( RNX ) 

• 

A2F31660 

EN 1 1 ( 32 I 

ENA (0) 

9 

A2F31670 

STAUWRB + 21 ) 
STA (WRB*18> 

I JR 1 1 L ) 

•CLEAR LAST OF BUFFER 

9 

A2F31680 

A2F31690 

CN 1 3 ( 11 > 

ENI 2 ( 95 ) 

.•LOOP 

A2F3L700 

EN I 1 ( 7 ) 


• * TO 

A2F31710 

IDA2 ( RDB ) 

LRS 16) 

.•PACK 

A2F31720 

I N I 2 ( -1 ) 

IJPHL-1) 

.♦THE 

A2F31730 

STQ3 1 WRB+ 1 ) 

l J P 3 ( L“ 3 ) 

.♦title 

A2F31740 

LDA ( NSS ) 

STA (WRB+13) 

9 

A2F31750 

LDA (1.) 

STA (WRB+14) 


A2F31760 

LDA <1. ) 

STA (WRB+15) 

9 

A2F31770 

LDA (NCHAN) 

STA (WRB+16) 

9 

A2F31780 

LDA (51.2) 

STA ( WRB+17 ) 

• 

A2F31790 

ENA ( 1 ) 

STA (WRB+19) 

• 

A2F310OO 

LDA ( SMOTHD ) 

STA (WRB+20) 

♦ 

A2F31810 

LOA (FIRST) 

STA ( WRB ) 


A2F3102O 

EN I 6(0) 

SL J4 ( OS EL ) 

•SELECT A TAPE 

A2F31830 

LDA (OTAR) 

SAU (L+l) 

• 

A2F31840 

EXF (*) 

EXF (11100B) 

•INTERRUPT ON END OF WRITE 

A2F31B50 

ENA (Tl) 

SAL (Z+l) 

9 

A2F31860 

EXF1 (TO) 

EXF ( 50 34 IB ) 

.START TYPEWRITER AND CONVERTER 

A2F 31870 

EXF5( Ai) 

* GET FIRST SCAN * 

EN 1 4 ( 0 ) 

.START CHANNEL 5 

A2F31880 

A2F31890 

A2F31900 

A2F31910 

A2F31920 

A2F31930 

EXF7 ( 108 1 

SLJ (OSTXJ 

.EXIT IF TYPEWRITER DONE 

A2F31940 

EXF7 ( 5032AB 1 

SLJ (L-l) 

• 

A2F31950 

L I U6 ( Z+5 ) 

EXF ( 50 34 1 B ) 

.GET ADDRESS, START CONVERTER 

A2F31960 

EXF5 (Al ) 


•RESTART CHANNEL 

A2F31970 

OAGN RAO (SCANS) 
STA (LOST) 

STA (BURNED) 

.COUNT SCAN 

* 

A2F31980 

A2F31990 

1AGN ENA6 ( * ) 

AJP3I L-5) 

.REJECT SCAN IF TOO SHORT 

A2F32000 

A2F32010 




A2F32020 

* INITIALIZE ADDRESSES * 


A2F32030 




A2F32040 

A2F32050 

1 1 A S IU6 ( 3TS T } 

ENA6 ( * ) 

* 

A2F32060 

SCM (-00) 

SAU (1IN5I 

.TEST Al UNDERSIZE 

A2F32070 

2 I A INA (•) 

SAU ( IT ST J 

.TEST Al OVERSIZE 

A2F32080 

1NA ( A 1+ 1 ) 

SAU ( 2 I N6 ) 

.DISPLACEMENT FDR STRETCHING A2 

A2F32090 

SCM (-00) 

SAU 11IN6) 

.DISPLACEMENT FOR STRETCHING Al 

A2F32100 

31 A INA (*) 

MU I (5) 


A2F32110 

DVi (NCHAN) 

STA ( HR 0 + 1 3 ) 

.NUMBER OF SAMPLES 

A2F32120 
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STA 

(NSS) 

MU I 

INCHAN) 

• 



A2F32 1 30 


CVI 

(5) 



• 



A2F32140 


CJP 

<L+1 > 

INA 

(1) 

.NUMBER OF 

WOROS PER LINE 


A2F32150 


STA 

[SIZE) 

SAL 

(5NXT) 

♦ 



A2F32160 


SAL 

( 1NXT ) 

INA 

(-1) 

* 



A2F32170 


SAU 

(UNX) 

SAL 

i UNX) 

* 



A2F32180 


SAL 

(51 A) 

INA 

(All 




A2F32190 


SAU 

( IIN7) 

INA 

(1) 

•LIMIT FOR 

STRETCHING Al 


A2F32200 


SCM 

(-00) 

SAU 

( 1IN4) 

• 



A2F32210 

41 A 

SAL 

f IK) 

ENA6 ( * ) 

.TEST FOR 

STRETCHING NEEDED IN Al 

A2F32220 


SAU 

( 2FLG ) 

SAU 

( 1IN1) 

•START OF 

A2 


A2F32230 


SAL 

( 1IN3) 

SAL 

( 2 I N 2 ) 

* 



A2F32240 

5 I A 

SAU 

( IK) 

INA 

1 * ) 

9 



A2F32250 


SAU 

( 2 I N7 ) 

INA 

( 1) 

.LIMIT FOR 

STRETCHING A2 


A2F32260 


SCM 

(-00) 

SAU 

(2IN4) 

.TEST FOR 

STRETCHING NEEDEO IN A2 

A2F32270 


ENA6 

(0) 

INA6I0) 

9 



A2F32280 


SCM 

t-OB) 

INA 

( Al) 

• 



A2F32290 


SAU 

(2IN5) 

SCM 

( -OB ) 

.TEST FOR 

A 2 UNDERSIZE 


A2F32300 

61 A 

INA 

(*) 

SAU 

( 1Q0M) 

.START OF 

M 


A2F32310 


SAU 

( 10 1 M ) 

SAU 

( 2 OOM ) 

• 



A2F32320 


SAU 

( 2 I Nl ) 

SAL 

( 2 [ N 3 ) 

* 



A2F3233Q 


SAU 

(OK) 

ADO 

(SIZE) 

* 



A2F32340 


SAU 

( LOON) 

SAU 

( 1 0 IN ) 

.START OF 

N 


A2F32350 


SAU 

(200) 

ADO 

(SIZE) 

• 



A2F32360 


SAU 

(7WRT) 

ADD 

(SIZE) 

•START OF 

TAPE CORRECTION AREA 

A2F32370 


SAU 

( 6NXT ) 

SAU 

(OWRT) 

•START OF 

FIRST OUTPUT BUFFER 


A2F32380 


SAL 

( 6WRT ) 

SAU 

(2NXT) 

* 



A 2F 32 390 


INA 

(-1) 

SAL 

l 3K ) 

•LAST WORD 

OF TAPE CORRECTION 

AREAA2F32400 


SAU 

(onxt) 



• 



A2F32410 


WAIT FOR GOOD SCAN 



EXF7 ( 10B ) 

SLJ I OS TX ) 


EXF7 1 50324B ) 

SLJ (L-l) 


LIU6IZ+5) 

EXF (50341B) 

1 TST 

ENA6 ( * I 

AJP2 I 2B IG > 

2T$T 

INA (•) 

A JP2 ( OK ) 


EXF5 ( Al ) 

SIU6U+3) 


ENA (0) 

S TA6 ( Z ) 

3TST 

I SK6 ( * ) 

SLJ IL-1) 


EN I 6 ( * ) 

SLJ ( 0 AGN ) 

2BIG 

£XF5(A1> 

SLJ ( 0 AGN ) 

* GOOD SCAN, START 

EVERYTHING • 

OK 

EN 1 5 ( * ) 

S IL 5 ( Z + 5 ) 

IK 

EXF5 ( * ) 

ENA6I *) 


AJP2IL+1 1 

RAO (STRETCH) 

2 K 

RAO (SCANS) 

ENA6( * ) 


AJP3U + 1) 

RAO (LONG) 

3K 

ENA (-0) 

STA [#) 


ENA (OINT) 

6XF4 ( WRB ) 


EXF (50 342B } 

EXF (10B) 


SAL (Z+70) 

SLJ I 1FLG ) 


.EXIT IF TYPEWRITER DONE 
.WAIT FOR END OF LINE 
.SAVE SIZE, RESTART CONVERTER 
.JUMP IF OVERSIZE 
.JUMP IF SIZE IS OK 


A2F32420 
A2F32430 
A2F3244Q 
A2F32450 
A2F32460 
A2F32470 
A2F32480 
A2F32490 
A 2F32500 
A2F32510 


» TOO SMALL, START OVER, SAVE SIZE A2F32520 


.CLEAR LEFTOVER DATA 
.LOOP I CLEAR AN EXTRA WORD) 
.RESTORE SIZE, GO RE- I NI T I AL I Z E 
.TOO BIG, GO RE- INITIALIZE 


A2F32530 
A2F32540 
A2F32550 
A2F32560 
A2F32570 
A2F32580 
A2F32590 
A2F32600 
A2F32610 
A2F32620 
A2F32630 

.IF SCAN WAS SHORT, COUNT STRETCH A2F32640 
.COUNT SCAN A2F32650 
.IF SCAN WAS LONG, COUNT A TRIM A2F32660 
.STORE THE TAPE CORRECTION WORD A2F32670 
.START THE TAPE A2F32680 
.SELECT CONVERTER AND CONSOLE INTSA2F3269Q 
.CONNECT THE INT. ROUTINE AND GO A2F32700 

A2F32710 


.START BUFFER FOR NEXT SCAN 
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PROCESS THE DATA I FROM ADTEST2 ) 
(INDENTED CARDS ARE UNMODIFIED) 
(MAY NCT USE INDEX 6) 

I SMOOTHS OF 3 THROUGH 16) 


GET INPUT BUFFER TC EMPTY 


LOO 

I SK7 ( OFLG ) 

SLJ (L+3) 


1 FLO 

ENA (Al) 

ISK {*> 



SAL (OFLG) 

SLJ l L + 3 ) 


2FLG 

ENA 1*1 

ISK (*> 


OF LG 

SIL (OFLG) 

ENI (*) 



EXF7 1 LOR 1 

slj (ooun 


1INX 

ENI1(*) 

ENI2l*l 


* FIRST PASS - STORE 

VALUES * 



I JP41 1 QOC ) 

S AU ( L + 1 ) 

.SET 


LDA2 ( * ) 

SCM(MASK ) 

. B I A 

100S 

ARS ( • ) 

LDQ (MASK2) 

ft 

100M 

ST L2 ( * ) 

LDO (MASK It 

ft 

100N 

ST L2 1 * ) 

1 JP2IL-3) 

ft 


SLJ (100D) 


ft 

* REMAINING PASSES - 

ADD TO STORAGE * 

100C 

I N 14 ( 1 ) 

S AU ( L+ 1 ) 

.SET 


LD A2 ( * ) 

SCM(MASK) 

.eiA 

LOIS 

AftS (*) 

STA (DATUM) 



LOG (DATUM) 

LDL (MASK2) 

• 

L01M 

RAD2 ( » ) 

LDL (MASK1) 

. 

LOIN 

R AD2 ( * ) 

IJP2IL-A) 

• 

* RELEASE THE INPUT 

BUFFER * 


100D 

ISK7 (OFLG) 

SLJ (L + 2) 

. 


S I LI ( 2F LG i 

SLJ ( L+ 2 ) 

• 


Sill (1FLGI 

ARS (0) 

. 

1ULE 

ISKM*) 

SLJ (100) 

•( 


WHICH BUFFER 
WAIT FDR BUFFER 1 

WAIT FOR BUFFER 2 

EXIT IF DONE 


A2F32720 
A2F32730 
A2F32740 
A2F32750 
A2F32760 
A2F32770 
A2F32780 
A2F32790 
A2F32800 
A2F328 1 0 
A2F32820 
A2F32830 
A2F32840 
A2F32850 
A2F3206D 
A 2F3207O 
A2F32880 
A2F32890 
A2F32900 
A 2F 32910 


CONTINUE 


* GET EMPTY OUTPUT BUFFER * 



ENA7I7NXT) 

INA 

( 1) 


5NXT 

SAU 

(200S) 

INA 

< *) 



SAL 

( 7NXT ) 

SAL 

( L + 

3) 


INA 

l -37777 R ) 

ajp: 

ML + 

2) 

6NXT 

ENA 

(*) 

SLJ 

(L- 

3) 


- 


I SK7 ( * ) 



PLACE BUF SIZE IN INDEX R EG I STERS A2F32920 

A2F32930 

A2F32940 

A2F32950 

A2F32960 

A2F32970 

UP A2F32980 

S DATA UPWARD BY 5.0 VOLTS A2F32990 

.SHIFT RIGHT FOR MAXIMUM CAPACITY A2F33000 

A2F33010 
A2F33020 
A2F33030 
A2F 33040 
A2F33050 
A2F33060 
A2F33070 
A2F33080 

UP A2F33090 

S DATA UPWARD BY 5.0 VOLTS A2F33100 

SHIFT RIGHT FOR MAXIMUM CAPACITY A2F33110 

BRING EVEN VALUES INTO ACCUMUL ATOA2F3 3 1 20 
ODD VALUES INTO ACCUMULATOR A2F33130 

CONTINUE A2F33L40 

A2F33150 
A2F33160 
A2F33170 
A2F331 80 
A2F33190 

WHICH BUFFER A2F33200 

RELEASE BUFFER 2 A2F33210 

RELEASE BUFFER 1 A2F33220 

CONTINUE A2F33230 

A2F33240 
A2F332 50 
A2F33260 
A2F33270 
A2F33280 
A2F33290 
A2F33300 
A2F33310 
A2F33320 
A2F33330 
A2F33340 
A2F33350 


.STORE START OF NEW BUFFER 
.STORE END OF NEW BUFFER 
.DOES NEW BUFFER EXIST 
.NO-GET FIRST ONE 
.WAIT FOR BUFFER EMPTY 
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066££dZV 
086EC3ZV 
0i6££3ZV 
096EE3ZV 
0£6EEdZV 
0^6££d2V 
0£6££dZV 
0Z6££dZV 
Q16£EdZV 
006E EdZV 
068££dZV 
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* BUFFER 1 * 

*#•*«#***•#» 



A2F34000 

A2F34010 

A2F34020 

A2F34030 

1 1 N 1 

ENA <*) 

SAL (Z + 5) 

.SET ADDRESS OF END OF BUFFER 

A2F34040 


I SK7 ( 1 FLG ) 

SLJ (L + 4) 

♦IS THE BUFFER EMPTY 

A 2F 34050 


EXF 1503416) 

SLJ4< 0A20) 

..FULL-START AND GD WAIT 

A2F3406Q 

11 NO 

KAO (LOST) 
SLJ (L-31 

RAO (SCANS) 

..COUNT LOST SCAN 
..GET NEXT SCAN 

A2F34070 

A2F340B0 


EXF ( 5034 1 B ) 
S L J4 { 0A2D ) 

EXF5I A 1 ) 

.EMPTY-START CONVERTER AND CHANNE LA2F34090 
.GO WAIT FOR END OF LINE A2F34LOO 


EXF7 ( 5 l B ) 

SLJ ( L + 2 ) 

.TEST FOR CHANNEL ACTIVE 

A2F34U0 


RAG (LONG) 

SLJ ( 1 I N ) 

.NOT ACTIVE-LONG SCAN-ACCEPT IT 

A2F34120 

1 IN3 

L l U6 ( 2 ♦ 5 ) 

EXF5 ( • > 

.SAVE ADDRESS-STOP CHANNEL 

A2F34130 

1 1 N4 

ENA6 l * ) 

AJP2I 1IN) 

.TEST FOR FULL SCAN 

A2F34X40 

1IN5 

ENA6 ( * ) 

A JP 2 ( L <■ 2 ) 

.SHORT. ..STRETCH OR SCRATCH 

A2F34150 


ft AO (SHORT) 

SLJ (1INQ) 

.♦SCRATCH IT 

A2F34160 

1IN6 

LDA61*) 

S TA 6 ( Z ) 

..STRETCH IT 

A2F34170 

1 1 NT 

I SK6 l » ) 

SLJ IL-l) 

. .CONTINUE 

A2F34180 


RAC tSTRETCH) 

ARS (0) 

• • 

A2F34190 

1IN ft AO (SCANS) 

* BUFFER 2 * 

o * «* *«*«»*«* 

SIL ( 1 FLG ) 

.ACCEPT THE SCAN-RELEASE THE BUFF 

A2F34200 

A2F34210 

A2F34220 

A2F34230 

A2F34240 

A2F34250 

2 I NX 

ENA (*) 

SAL (Z + 5) 

♦SET ADDRESS OF END OF BUFFER 

A2F34260 


ISK7I2FLG) 

SLJ (L+4) 

.IS THE BUFFER EMPTY 

A 2F 34270 


EXF ( 50 34 IB ) 

SLJ4(0A2D) 

. .FULL-START AND GO WAIT 

A2F34280 

2 1 NO 

RAO (LOST) 
SLJ (L-3) 

RAO (SCANS) 

..COUNT LOST SCAN 
..GET NEXT SCAN 

A2F34290 
A 2 F 34 300 

2IN2 

EXF ( 50 34 1 B ) 
$ L J4 ( 0A20 ) 

EXF5( *> 

•EMPTY-START CONVERTER AND CHANNEL A2F34310 
.GO WAIT FOR ENO OF LINE A2F34320 

OIN 

EXF7 ( SIB ) 

SLJ (L + 2) 

.TEST FOR CHANNEL ACTIVE 

A2F34330 

2IN3 

RAO (LONG) 

SLJ ( 2 I N ) 

-NOT ACTIVE-LONG SCAN-ACCEPT IT 

A2F34340 

LIU6IZ+5) 

EXF 5 ( * ) 

.SAVE ADDRESS-STOP CHANNEL 

A2F34350 

21N4 

E N A 6 ( * ) 

AJP2( 2IN) 

.TEST FOR FULL SCAN 

A2F34360 

2 I N5 

ENA6 ( *) 

A JP 2 ( L + 2 ) 

.SHORT.. .STRETCH OR SCRATCH 

A2F34370 

2 I N'6 

HAD (SHORT) 

SLJ ( 2 I NO ) 

..SCRATCH IT 

A2F34380 

LDA6 ( * } 

S T A 6 ( Z } 

..STRETCH IT 

A2F34390 

21N7 

ISK6I » J 

SLJ (L-l) 

..CONTINUE 

A2F34400 


RAO (STRETCH) 

SLJ ( 2 t N ) 

• • 

A2F34410 

A2F34420 

* 


t« 


A2F34430 

A2F3444Q 

* TAPE 

• 

: WRITING LOGIC 

« 

* 


A2F34450 

A2F34460 



Ml 


A2F34470 

A2F34480 

A2F34490 

* GET 

-**#S'**#*#**#«#-S 



A2F34500 

A BUFFER OF OUTPUT * 


A2F34510 

ONXT 

ENA (*) 



A2F34520 

A2F34530 

INA (1) 

♦ 

A2F34540 

1 NXT 

SAU (1WRT) 

INA (■) 

.STORE START OF NEXT BUFFER 

A 2F 34550 


SAU ( 9NXT ) 

SAU (OWRT) 

•STORE END DF NEXT BUFFER 

A2F34560 


SAU (ONXT) 

SAU l 4NXT ) 

. 

A2F34570 

2nxt 

INA I -37777 B ) 

A JP3 ( 3NXT ) 

.DOES NEW RUFFER EXIST 

A2F34580 

ENA ( * ) 

SLJ ( 1NXT ) 

.ND, CHANGE TO FIRST ONE 

A2F34590 

3 NXT 

ENA ( 20B ) 

SIU (8NXT) 

.CLEAR FORCED INTERRUPT 

A2F34600 

4NXT 

I SK7 ( * ) 

SLJ (OWRT) 

.IF NEW BUFFER FULL, GO TO WRITE 

A2F34610 


SAU (8NXT) 
SLJ ONXT) 

SLJA(OTAP) 

.WAIT FOR BUFFER FULL INTERRUPT 
•GO TEST IT AGAIN 
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* WRITE A RECORD * 


OWRT 
1 WRT 
2 HR T 
3WRT 
<tWRT 
9NXT 


ENA (*> 
E XF4 ( * ) 
EXF7 < * ) 
EXF7 ( * ) 
EX F7 ( * ) 
SIL 1*1 


SAL (Z + 4) 
SL J4( or AP ) 
SLJ (5WRT) 
SLJ ( SWRT ) 
SLJ4(OEOT) 
SLJ ( ONXT ) 


• CORRECT TAPE ERRORS 


SWRT 

EXF7 { * ) 

SLJ 

( 9WRT ) 

6 WRT 

EXF (•) 

ENA 



SAL ( Z-t-4 ) 

SL J4 ( OBKS ( 

7WRT 

EXF4 ( • ) 

SLJ4 ( QT AP) 

8 WRT 

EXF? ( * J 

SLJ 

( 9WRT ) 

8WRT1 

EXF (*| 

SLJ4( OBKS) 


RAC (ERRORS) 

SLJ 

(OWRT J 

9 WRT 

SLJ4IOBKS) 




SLJ4 (OEOT ) 




RAC (ERRORS) 

SLJ 

(OWRT) 


* TAPE CONTROL SUBROUTINES 


* BACKSPACE ONE RECORD • 


1 OKS 
0S5KS 


EXF 

SLJ 


(*} 

1*1 


SLJ4 t OTAP ) 
SLJ (L-l) 


.SET END 

♦START WRITE, WAIT FOR DONE 
.TEST FOR PARITY ERROR 
.TEST FOR LENGTH ERROR 
.CHANGE TAPE IF END OF REEL 
•RELEASE BUFFER, GO GET NEXT ONE 


•TEST FOR END OF TAPE 
•SELECT BCD, SET END OF CLEAR BUFFA2F348I0 


A2F3465Q 

A2F34660 

A2F34670 

A2F34680 

A2F34690 

A2F34700 

A2F34T10 

A2F34720 

A2F34730 

A2F34740 

A2F34750 

A2F34760 

A2F34770 

A2F34780 

A2F34790 

A2F346O0 


.BACKSPACE 
.WRITE BLANK TAPE 
.TEST FOR END OF TAPE 
.NOT END, BIN MODE AND BACKSPACE 
.COUNT ERROR, GO RE-WRITE 
.BACKSPACE 
.CHANGE TAPE 

.COUNT ERROR, CO RE-WRITE 


.BACKSPACE/WA I T READY 
.EXIT/ENTRY 


* END OF TAPE - UNLOAD TAPE AND GET A NEW ONE 


OEOT 

SLJ <*l 

EXF (*) 


ENA (SPACES) 

SAL (Z+4) 


EXF4 ( ENDREEL ) 

SLJ4I0TAP) 

1 EOT 

EXF (•) 

EN I 6 ( * ) 

2EOT 

1SK6(*> 

ARS (0) 

3 EOT 

SI L6 ( 1 EOF ) 
SLJ IOEOT) 

SLJ410SEL) 


.EXIT/ENTRY, SELECT BCD 
.SET END OF EOT RECORD 
.WRITE ENDTAPE RECORD/WAIT READY 


A2F34820 

A2F34830 

A2F34840 

A2F348 50 

A2F34860 

A2F34870 

A2F34880 

A2F34890 

A2F34900 

A2F34910 

A2F34920 

A2F34930 

A2F34940 

A2F34950 

A2F34960 

A2F34970 

A2F34980 

A2F34990 

A2F35000 

A2F35010 

A2F35020 

A2F35030 

A2F35040 

A2F35050 

A2F35060 

A2F35070 

A2F35080 

A2F35090 

A2F35100 

A2F35U0 


.UNLOAD TAPE, PUT COUNTER IN I NDEX A2 F35 1 20 


» SELECT TAPES * 


9SEL 

BS5 

(8) 


OSEL 

SLJ 

<*) 

LDA6 ( 9S EL ) 


SAU 

( L + l ) 

SCL (77B) 


EXF 


SAU (1TAP) 


INA 

m 

SAU (5SEL) 


SAU 

( 8WRT1) 

INA (I) 


SAL 

(OEOT) 

SAU (6WRT) 


.CYCLE TAPE COUNTER 
•STORE COUNTER, SELECT NEW UNIT 
•RETURN 


.TAPE SELECT CODE TABLE 

.EX I T /ENTRY , GET SELECT CODE 

.STORE XCOU1 

.STORE XCOOO, SELECT TAPE 
•STORE XC001 


.STORE XC002 


A2F35I30 

A2F35140 

A2F35150 

A2F35160 

A2F35170 

A2F351B0 

A2F35190 

A2F35200 

A2F35210 

A2F3S220 

A2F35230 

A2F35240 

A2F35250 

A2F35260 

A2F35270 
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INA 

( 1) 

SAU 

( 2WRT) 

•STORE XC003 

A2F3520O 


INA 

(1) 

SAL 

IOT API 

•STORE XC004 

A2F3529D 


INA 

m 

SAU 

t 3WRT) 

•STORE XC005 

A2F35300 


INA 

( n 

SAU 

( 1BKS 1 

•STORE XC006 

A2F35310 


I NA 

m 

SAU 

t ieot) 

•STORE XC007 

A2F35320 


SAU 

( 4WR T ) 

SAU 

( 5WRT ) 

m 

A2F3S330 


SAU 

( 8WRT I 

ARS 

(0) 

• 

A2F35340 

5SEL 

EX F7 ( * } 

SLJ 

(OSEL) 

.RETURN IE TAPE READY 

A2F35350 







A2F35360 







A2F35370 

* 


« 




A2F35380 

* END 

THE 

RUN * 




A2F35390 

* 


* 




A2F35400 


NEW TAPE NOT RE AC Y , CRASH THE RUN * 


9EGT 


EXF (11B) 

EXF7 c 110 I 
LDA ( Z + 5 ) 

EXF5 [*) 

SAU IZ+7) 

ENI (0) 

ENA (WCNKITY) 
ENA (*) 
SLJ4(0Z) 

SLJ ( 1 £ NT ) 


EXF (500000) 
EXFUT1 ) 

SAU (L+l) 

ENA (L+2) 

SLJ ( Z+7 ) 
ENKO) 

SL J 4 ( MT R ) 

SAL (Z+7) 


AFTER INCOMPLETE START # 


OSTX EXF5 ( * ) 
SLJ4I0Z) 
SLJ ( 1ENT ) 


EXF (50000B) 


AFTER NORMAL RUN * 


OOUT 

LIUK8NXT) 

IJPliL+2) 


SLJ 

(L“l) 



LDA 

(OTAP) 

SCL (770) 


SAU 

(L+ii 

SAL (L+l) 


EXF 

(*) 

EXF7 ( * ) 


INA 

(3) 

SAU ( 20UT ) 


INA 

(4) 

SAL ( 30UT ) 

10UT 

ENA 

(*) 

SAL (Z+7) 

2 OUT 

EXF 

<*) 

SLJ A ( OZ ) 


ENA 

(DONE) 

SLJ4I MTR) 


ENA 

(9) 

SLJ41RNX) 


LDA 

( RDB ) 

INA (-300) 

30UT 

AJP1 ( 1 ENT ) 

EXF (*) 


SLJ 

( 1 ENT ) 


* PRIM STATISTICS 

• • 

* 

• * 


OZ SLJ (*) 

EN 1 5 ( NUMS T I X ) 


CLEAR INTERRUPT AND CONVERTER 

STOP CHANNEL 5 

LEAVE INTERRUPT MODE 
00 NOTHING TO CATCH JUMP 
WRITE * NO TAPES* MESSAGE 

PRINT STATISTICS 
EXIT 


•CLEAR CHANNEL 5 
•PRINT STATISTICS 
• EXIT 


WAIT FOR LAST TAPE WRITE 


CLEAR TAPE INTERRUPT SELECT 

STORE EOF SELECT CODE 

STORE UNLOAD SELECT CODES 

RESTORE INTERRUPT 

WRITE EOF i GO PRINT STATISTICS 

ASK IF USER WANTS TAPE UNLOADED 

GO GET REPLY 


UNLOAD TAPE IF YES 


•EXIT/ENTRY 


ENA (0) 


A2F35410 
A2F35420 
A2F35430 
A2F35440 
A2F35450 
A2F35460 
A2F35470 
A2F35480 
A2F35490 
A 2F 35 500 
A2F35510 
A2F35520 
A 2F 355 30 
A2F35540 
A 2F35550 
A2F35560 
A2F35570 
A2F35580 
A2F35590 
A2F35600 
A2F35610 
A2F35620 
A2F35630 
A2F35640 
A2F35650 
A2F35660 
A2F35670 
A2F35680 
A2F35690 
A2F35700 
A2F35710 
A2F35720 
A2F35730 
A2F35740 
A2F35750 
A2F35760 
A2F35770 
A2F35780 
A2F35790 
A2F350OO 
A2F35810 
A2F35820 
A2F35830 
A2F3504O 
A2F35850 
A2F35860 
A2F35870 
A2F35880 
A2F35890 
A2F35900 
A2F35910 
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z 

THSSICRRCRS ) 

SLJ (OZ) 

.SEARCH/EXIT 

A2F35920 


ENA5 { J J ) 

I NA5 ( 0 ) 

.GENERATE ADDRESS 

A2F35930 


I NA5 i 0 ) 

LDQ5( ERRORS ) 

.GET VARIABLE 

A2F35940 


SAL (3Z) 

SAL (4Z) 

•SAVE ADDRESS 

A2F35950 


GJP3I5Z > 

LDA (SPACES) 

.JUMP IF MINUS, GET BLANKS IF 

NOT A2F35960 

21 

ST A7 ( 4Z ) 

E N I 2 ( 6 ) 

.SET BLANKS 

A2F35970 

L D A5 ( ERRORS > 

ENQ (0) 

. 

A2F35980 


DVI (10) 

S T A 5 ( ERRORS > 

.GET LEAST SIGNIFICANT DIGIT 

A2F35990 


GJPKL+l) 

ENQ (12B) 

.SET TO 12B IF ZERO 

A2F36000 

3Z 

IDL ( 77B ) 

LDQ <*> 

.PACK DIGIT 

A2F36010 

4Z 

LRS (61 

STU («) 

• 

A2F36020 


LDA5 < ERRORS ) 

AJP (L+2) 

•EXIT IF NUMBER DONE 

A2F36030 


I JP2(2Z) 

SLJ ( 5Z ) 

.LOOP/SET MANY IF OVERFLOW 

A2F36040 


LDA (SPACES) 

LLS (6) 

* 

A2F36050 


LRS (6) 

I JP2 ( L ) 

•FILL WITH LEADING BLANKS 

A2F36060 


STQ7I4Z ) 

SLJ { L ♦ 2 ) 

.STORE AND GO PRINT 

A2F36070 

5Z 

LDA ( MES6 ) 

STA7( 4Z ) 

.SET / MANY / 

A2F36080 


ENA7KZ ) 
SLJ (1Z) 

END 

S L J 4 ( M T R ) 

•PRINT LINE 
.LOOP FOR MORE 

A2F36090 

A2F36100 

A2F36H0 

A2F36120 


PROGRAM CSClt I ARG1 » I ARC 2 . IARG3) 

C NEEDS PROCESS, UNPACK3, 8UFPACK. 

C INPUT ON 3, OUTPUT ON 4. 

C PROGRAM HAS 'CONTINUE FILE* OPTION 
C ‘TIME* CHANNEL SHOULD BE LAST CHANNEL ON TAPE. 

C DESIGNEC PRIMARILY FOR DOING DYNAMIC CLAMPING TO THE DARK LEVEL, 

C AND DYNAMIC SCALING TO THE SUN SENSOR. 

C SCALE FACTORS ARE COMPUTED ONLY WHEN AIRCRAFT ROLL (AS INDICATED 
C 8Y THE SIGNAL IN A DIGITIZED 'TIME* CHANNEL) IS WITHIN CERTAIN 
C LIMITS. LIMITS WILL BE COMPUTED BY THE PROGRAM. 

C HAS ’DESKEWING* CAPABILITY (’SKEWS = • 0 FOR NO DESKEWING) 

C THE 'CLAMPING REGION’ WILL NOT BE DESKEWED PRIOR TO CLAMPING DATA 
C ’SCALING REGION’ IS DESKEWED PRIOR TO SCALING (MINIMUM SKEW .EQ. 0) 
C 

C SAMPLE REGION LINE NO.* 0 TO SKIP COMPUTATION OF LIMITS 
C A SCALE FACTOR OF * 0. • WILL SIGNAL THAT PREVIOUS VALUE SHOULD BE 
C USEO (* 1.0 » WILL BE USED IF NO PREVIOUS VALUE WAS ENTERED). 

C ONLY ONE * 0.’ NEEDS TO BE ENTERED IF ALL SCALE FACTORS ARE TO 

C BE LEFT UNCHANGEC. 

C NO. OF FILES TO SKIP = -0 TO TERMINATE THE RUN 

DIMENSION ISKS 1 12) 

DIMENSION I VSSt 12 ) 

DIMENSION l KEY ( 6 ) , KE Y WD ( 6 ) 

DIMENSION D A T A ( 6000 ) , $HIFT(5) 

DATA ( KEYWDI 0 ) =00004546 2 37 14465B , C04 546475 1 7 1 4523B , 2020202020202020 
1 B, 204546634 36 14447B, 204 546 2 263614 365B, 0000000 0006 346458) 

COMMON I CODE ( 12 ) , ISKEW113), 1VCI12), IVSI12), NHK12), NL0(12), 

1 USC ( 12 ) , VC ( 1 2 ) , VCL ( 12 ) , VSCI12), VSSI12), WSSI12) 

COMMON DEL (51, ID(12), JVI50), P0SI5), SAMP(500I, X(50) 

CAT A ( SH I FT = 9, 18, 27, 36, 45), (MAX = 777BI, (MIN = -400B ) 

COMMON V 1 6000 ) , RESERVEI12), IDl, ID2, BANG, DANG 

COMMON L90, KEY, NPTS, NUNES, NSA, NSB, KS, NA, N6, KP, NEWR, 

1 TITLEI12), IAGK12), TAG2I12I, NSS, NCHAN , KR, CDNV, IPOS, 

2 I PACK, INT 

EQUIVALENCE ( DATA ( 1 ) , V(l)), (MINSKEW, ISKEW(O)) 

INTEGER DATA, DEL, POS, SAMP, SHIFT, TAG 1 , TAG2, TITLE, VC, X 
LOCI R LB = 13, TEST = 70, ZRO = 0) 

PARTMAP 
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CALL CCRECON 
N5 = 5 

CALL ZERO I V SC , VSCI 12), 1.) 

I CRASH = 1 
FF = 59.0 

I END = 0000000000654564B 

LDAIIARG1) SUB(IENO) AJP0I44) . JUMP TO END RDUT1NE 

LIU11TESTI I JP1 ( 1 ) SLJI2) 

1 WR I T E ( 1 9 , 10 1 ) 

101 FORMAT!/ * EXECUTION BEGINS * /) 

SLJIL+2) . SKIP NEXT INSTRUCTION 

2 EXF(IOOOB) SLJ3I3) .IF NOT BATCH, START REAL-TIME CLOCK 

+ ENA I 77776 6 ) STAI IKEY) 

STAIIKEY+l) STAI IKEY + 5) 

ENA ( 00000 ) STAIIKEY+2) 

DO 500 1*0,2 

+ ENI6I 0C006B ) LDAl(IARGl) 

EQS6IKEYWD) SLJ(503A) 

ENA(OOOOi) STA611KEY) 

SLJI500) ENIOIO) 

503 A AJPOI5O0) STAIIKEY4-2) 

500 CONTINUE 

ENA 1 00000 ) STA(IARGl) 

STAIIARG2) STAIIARG3) . ZERO OUT PAR FIELDS 
+LDAIIKEY+S) AJP3I3) 

ENA ( 00000 ) STAI IKEY+2) 

3 CALL PROCESSIDATA, 1) 

LDAIKEY) SUB ( 5 1 AJPI91) 

CALL ZERO I W SS , WSSI 12) ) 

ENAI77776B) STAI IKEY+3 ) 

STA I IKEY+4 ) ST A t I KEY +6 ) 

IF IRST=1 

Nk I T E ! 9 , 107 I 37 B 

107 FORMAT I • NEW TITLE* R 1/ ) 

LDAINCHAN) INAI77776B) 

+ SSKIIKEY) I N A 1 0000 1 ) 

STA(NC) 

IFINC .LT. 13)11,97 
11 CONTINUE 

READC9, 108 ) ID 

108 FORMAT < 12A8 ) 

LIUKTESTI I JP 1 1 2 1 1 SLJI22) 

21 WRITEI 9, 10B ) ID 

22 IFIIDI1) .EQ. 3HYES ) 11, 12 

12 1FIICI1) .NE. 2HN0 ) 13, 17 

13 DO 14 1=1,12 

14 TITLE! I ) = 10U) 

17 INT = 1 

IUNIT = 3 
IP0S2 = 1 
ISTART = 0 
LS TART = 0 
NOTITLE = 0 
KBAD = 0 
NB AD = 777B 

CALL ZERCINHI, ML Ot 12)) 

WRITE 19, 117) NC 

117 FORMAT!* SK EWS I • 1 2* I 5 ) = *) 

READ (9,103) ISKEW 
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103 FORMAT ( 151 5 ) 

L 1 0 1 t TEST ) IJPK17A) SLJ117C) 

17A WR I T E ( 9 f 1 03 ) USKEW(I), 1= 1 r NC ) 

17C MINSKEW = 10 
MAXSKEW = -10 


6 


105 

106 

1C9 


510 

520 

225 

511 


512 

514 


517 

513 

24 


CQ 6 J = 1 * NC 

MINSKEW = M I NO < M I N SK EW * ISKEW(J)) 
MAXSKEW - M AX 0 ( MAX SKEW* ISKEWU)) 
MAXSK E ini = MAXSKEW - M INSKEW 
WRI TE(9,105) 


F0RMATUCLAMP REGION# VCLAMP. 2 I 4 , 12F6. 3* / ) 

RE AC I 9 , 106 ) N ACL , NBCL , VCL 
FORMAT (214# 12F6.4) 

WRITE<9,109) 

FORMAT I * SC ALE REGION, VSCALE. 2I4,12F6 # 4 */) 
KE AC { 9 , 106 ) NASS# NBSS, USC 
NBCL- M I NO ( NBCL # NSS-MAXSKEW) 


NHSS = MIN0IN6SS, NSS-MAXSKEW) 
IFtNASS *G T * 7000 ) 510,511 
WRI TE ( 9 » 52 0 ) 

FCRMAT(*SMQCTHING CONSTANT = *) 
READ ( 9 T 225 ) FF 
FORMAT { F10 • 0 ) 

NASS = NASS-70 00 


LDA (NBSS ) 
AJPK512) 

SL J { 514 ) * 

SUB < NASS ) 
LOA(NBCL) 
IP0S2 = IPOS 
AJPH517) 
KBAO =-255 


ENCtOOOOlB ) 
STC ( IKEY + 4 ) 

A JP 3 (97) 
AJPLI 513) 


NBAD = 255 

STQI IKEY+3 ) SLJ( 24) 

SUB(NACL) AJP3I 97 ) 

NC ST = ( N ACL - 1 ) * NCH AN 
C0NV2 = 0.0 
FVS = 1.0/CCNV 
FFCCNV = FF + CONV 

DCONV 2 = 1.0 / ((1*0 + FF) * CONV) 

IKEY ( 6 ) = I KEY ( 4 ) + IKEY(5) - 1 
N = NC + 1 
LOA(O) E NC ( 00001 ) 

+THS6IISKEW) ENQ{ 77776B ) * CHECK FOR ANY DESKEWING 

STQ ( JSKEW ) 


C OETERMINE ACCEPTABLE LIMITS FOR SIGNALS [N 'TIME* CHANNEL 
NEWR = 0 

+ SSK ( IKEY+4) SLJ( 18) 

SSK(IKEY) SLJI18) 

LD A ( IKEY+2 ) A JP 0 ( IB ) 

STA(NSA) 

NSB= NSA+120 
I F I RST= 1 
KS= 2 

NA= NSS/3 

NB= NA + M I NO ( 50, NA) 

K P= I 

CALL ZERCIWSS, WSS( 12) ) 

194 FSS = NB - NA + 1 
NL - NSb - NS A 


57 



FORMERLY WILLOW RUN LABORATORIES, THE UNIVERSITY OF MICHIGAN 


2pl 


N L - NL/KS ♦ 1 
NP = 500/NL 

LIUl(TEST) IJP1136) SLJ(37) 

36 WRITE<9,103> N$ A , NSB, KS, N A , N8, KP 

37 ASSIGN 99 TO N90 
CALL SET EOF ( N90 J 

CALL FSKIPiOt NS A - NEWR » (UNIT) 

NEWR = NS A 
KM = KS - 1 

LDA( ZRC) SCL ( 7 777777777 77 7000R ) 

STA( ISRAN) 

CALL RAN1FI I SR AN ) 

SL J4 ( 320) .SET UP UNPACKING 

ISP = 0 
KR * 0 

00 94 NSCAN = N S A , NSB, KS 

DO 93 K - 1, KR 
NEWR = NEWR + 1 

LDAUUNIT) SLJ41RLB). READ ONE BINARY RECORD 

93 CONTINUE 

KR = KM 

NEWR = NEWR + 1 

SL J4 ( 300 ) . UNPACK CHANNEL NCHAN INTO X 

DO 94 J = l f NP 

K = I F I X ( RAN IF 1 - L ) * FSS - .0000011 + NA 
ISP * ISP + 1 

94 S A M P ( ISP) * X { K ) 

SLJI99A) 

99 NSB = NEWR 

99 A CALL SORT 1 ( SAMP, ISP, -1) 

NM ID = ISP/2 
NQTR = ISP/8 
NLCW = NM I C - NQTR 
NHIG = NM I D + NQTR 
MIN - SAMP(NLCW) 

MAX * SAMP(NHIG) 

41 ANGLE = 3 » ( MAX - MIN + 1) 

ANGLE = ANGLE/20. 

WRITEI9, 115) MIN, MAX, ANGLE 

115 FORMAT ( * INTEGER RANGE I * 13 * TO * 13 *) COVERS* F4.1 * DEGREES*/) 

42 CALL RST EOF 

18 CALL PR0CESS(DATA,2) 

NB = MINOINB, NSS-MAXSKEW) 

NB1 = MAXOINB, NBSS) 

NKP = KP* NCHAN - NC + { I SKEW ( 1 ) - I SKEW I NC ) ) *NC H AN 
LIUl(TEST) IJPK26) SLJI27) 

26 WR I T E ( 9 , 103 ) NS A , NSB, KS, NA, N8, KP 

27 NST = (NA - 1 * ISKEh(l) - MINSKEW) * NCHAN 
NPSS = (NBSS - NASSJ/KP + 1 

LDA( NOTITLE) AJP 1 ( 46 ) 

FP = KP 

FANG = < N A - 1) 

BANG = BANG + FANG*DANG 
GANG = FP *D ANG 
KSS = ( NB-NA ) /KP ♦ 1 
MSS = KSS*NC 

WRITEI4) TITLE, KSS, BANG, DANG, NC, CONV, IPDS2, 1, TAG1 , TAG2 
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♦ SSKUKEY+i) SLJI32) 

46 WRITEIN5, 114) TITLE. KSS, 8ANG. DANG. NC, CONV, IP0S2, NASS, 

I NBSS, NACL, NBCL, KP 

114 FORMATUH1 1 2 A8 / IB * PT S . /SCANL I NE , *7X*BANG= * F6 . 3 , 7X *DANG= * 

1 F6.3/ 110* CHANNELS, *7X*C0NV =*F4.l* I QUANT A / VOL T ) , * 7 X* I P0$ =* 

2 Il/8X*$TABLE SOURCE LOCATED BETWEEN POINT NUM0ERS*I5* AND *15 

3 *, CLAMPING REGION BETWEEN* 15 • AND * 15 *, EVERY* 13 /) 

32 VCIO) = 0 

WRITE19.108) 8H0PTI0NS .BHUSED 

DO 531 J = 0 , 5 

LDA2 ( I KEY ) AJP3IS31) 

WR1TEI9.108) KEYWO(J) 

531 CONTINUE 

L DA I I KE Y + 2 ) AJP1I532J 
Kft l TE ( 9 , 1 08 I 8HNO TIME , 8HSAMPL ING 
532 DC 4 J= 1 ,NC 

I CODE I J ) = J 

ISKEWINC-J+1) = ISKEWINC-J+l) - ISKEWINC-J) 

VC < J > = VCL I J ) *CONV + 0.5 
LDA2IUSC) A JPO I 4 A J 
VSC(J) = USC ( J ) 

GO TO 401 

4 A USC ( J ) = VSC I J ) 

401 VSStJ) * USCIJ) 

4 IVSSIJI = I VSSI J J )*FLOATI 2**15 J 



1 SKEW ( 1 ) = 0 



+ SSK ( I KEY* 1 ) 

SLJI 533) 


♦ SSKI IKEY+3 ) 

SLJI29) 


WRITE(N5,121) 

I VCL I J I , J= 1 , NC ) 

29 

SSKI IKEY+4) 

SLJI29Y) 


WRITE! NS, 121) 

(VSCIJ), J= l, NC 1 

2 9 Y 

SSKI IKEY+5 1 

SLJI29A) . 


SLJI533) 


29A 

WRITEIN5, 121) 

IVSSIJJ fJ~l,NC) 

121 

FCRMAT ( 7X, 12F6 

.3) 

533 

DO 7 J* 1 , NC 



ISKEW(J) = I SKEW ( J ) *NCHAN + 1 
7 ISKS(J) = I SKEW! J ) 

+SSKIICRASH) SLJI530A) 

CALL EXIT 
530A CONTINUE 

+SSM IKEY+3 ) SLJI541) . JUMP IF NO CLAMP 
NKPS = o 
NSTS*0 
NB1S = KSS 
NAS = 1 
KPS= 1 

NASS = INASS-NA+1) 

NSST = ( NAS S-l ) * NC 

CALL ZERO! ISKS, ISKSI 121,1) 

GO TO 542 

541 NSTS = NS T 
NKPS = NKP 

NB1S = N01 

NSST =(NASS-l)*NCHAN 
NAS = NA 
KPS = KP 

542 CONTINUE 

CALL BUFFON 
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CO 50 NSCAN = NSA,NSB,KS 
CALL PROCESS ( DAT A » 3) 

LOA {KEY) SU9I4) A J P ( 90 ) 

CALL ZERO! IVC, IVC< 12 ) ) 

JSC L = -1 

♦ SSKI l KEY+4 ) SLJ( 67) 

LDAI1KEY+2) AJP0I67) 

SSK ( I KEY ) S L J ( 6 7 ) 

N = INA-1 ) *NCF AN 

DO 60 1-1, KSS 
M = N 

N = M + NCHAN 
I DAT = DATA(M) 

JOAT a DATAIN) 

LDA(IDAT) AJPI60J 

SUfi(JDAT) AJP ( 65 ) 

60 CONTINUE 
GO TC 66 

65 I F I I CAT .LE. MAX .AND. IDAT .GE. MIN) 67, 68 
68 JSCL =1 

♦ SSKI IFIRST) SLJ( 50) 

67 K = 0 

♦SSK ( IKEY+3) SLJl 501 ) 

N = NCST 

DO 10 I P -NACL , N BCL 

DO 9 J=1,NC 
N - N + 1 

9 IVCl J) = I VC C J) + DAT A ( N ) 

♦SSKI IKEY J SLJl 10) 

N = N ♦ 1 
10 CONTINUE 

NVC = I NBC L - NACL ♦ 1) 

CD 30 J= 1 , NC 

30 IVCIJ) = IVCIJJ/NVC - VC ! J ) 


C CLAP? THE DATA 

+ SSK( IFIRST) SLJ4I 2STCL) 
N - NS T 

CO 35 I S = NA « NBl « KP 


DO 34 J- 1 , NC 
K = K ♦ 1 
N - N ♦ I SKEW I J > 


2CL 

LDA6IDATA) 

STA3(DATA) 

SUB2I IVC) 

.CLAMP NORMAL POLARITY 

DATA 

34 

+ SSK ( IKEY+6) 
CONTINUE 

SLJ4I 7C$T). 

IF CLAMP ONLY, COLLECT 

STATISTICS 

35 

N = N + NKP 




501 

♦SSKI IKEY+5) 

SLJl 16) 
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521 +SSK{ I K E Y + 4 } $LJ(502) 

SSK(JSCL) SL J ( 16 ) 

C COMPUTE SCALE FACTORS 

CALL BUFWAIT 

JVMAX = 0 

CALL ZERO! JV, JV( 50) ) 

N = NS ST 

CD 19 IP = 1, NPSS 

CO 19 J = 1, NC 
N ■= N + 1 

19 JVUP) = JV( IP) + OATA(N) 


DO 20 IP = 1* NPSS 
20 JVMAX = MAXOI JVMAX * JVC IP)) 

K = NPSS + 1 
LDAI JVMAX ) 

♦ EQS3UV) SLJI50J- LOCATE PEAK OF SUN SENSOR 

N * I NASS + <K-2)*KP)*NC 

DO 25 J = 1 f NC 
N = N + 1 

KVS = M SS ( J ) *CONV 2 
KSS ( J ) = KVS + DATA(N) 

WSS(J) = WSSIJMFVS 
AJPI25A) SL J 1250 

25 A VSSIJ) = 1, 

SLJI25) 

25C VSS(J) = VSCt J)/WSSt J) 

25 CONTINUE 

DO 347 1*1, NC 

347 IVSS(J) =* VSSIJ) * FLOAT! 2**15) 

C SCALE THE DATA 

16 K * 0 
NN “ NS T S 

DO 338 IS=NAS, NBlStKPS 

DO 340 J* 1 f NC 
NN = NN ^ ISKS(J) 

N * NN 
K * K + 1 

LCA6 C DATA) MUI2I IVSS ). 

ARS (15) ST A3 ( DAT A ) • SHIFT RIGHT AND STORE 

SLJ4{7CST) . 

340 CONTINUE 
338 NN * NN + NKPS 

505 IFIRST = -1 

C0NV2 = FFCCNV 
FVS = DC0NV2 
MMSCAN = NSC AN 
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+ SSKUKEY+1) SLJ1506). JUMP IF NOPRINT IS SPECIFIED 
SSMIKEY+3) SLJII5) . PRINT CLAMP VALUES IF CLAMPING 

WRITE(N5, 10 2 ) NSCAN, ( I VC ( J 1 * J=1,NC) 

102 FORMAT ( 1 H , 1316) 

MMSCAN = 0 

15 + SSMIKEY + 4) SL J ( 506 ) • PRINT SCALE FACTORS IF DYNAMIC SCALING 

WK I TE ( N5 j 1 19 ) MMSCAN, (VSS(J), J=1,NC) 

119 FORMAT I LH , I 5 , IX , 1 2F6 * 3 ) 

506 CONTINUE 

CALL BUFPACK(DATA,MSS) 

50 LDAIKEY) SUB l 3 ) AJPI90) 


90 WRITEI9, 125) 370 

125 FORMAT ( * CONTINUE FILE* RL> 

READ! 9, 108) IREP 
I F ( 1 REP «EQ. 3HYESJ43, 44 

43 NOTITLE * 1 
GO TC 18 

44 NOTITLE * 0 
ENOFILE 4 

+ SSKI IKEY+1) SL J t 3 > 

WR I T E ( N5 1 1 1 0 ) TITLE 

110 FORMAT ! 1H1 , 12A8) 

WR I T E ( N5 » 1 1 1 ) t I CODE ( J ) , J^l.NC) 

111 FORMAT ( 1H0 * CHANNEL* 1216/) 

WR1TE!N5,122) 8H * GT • , NBAC , I NH II J I ♦ J = 1 , NC ) 

122 FORMAT { A8 , 14,1216) 

WR I TE ! N5 , 1 22 ) 8H ♦ LT . , KBAD, ( NLO l J ) , J =1 , NC ) 

GO TO 3 

91 ENDFILE 4 
REWIND 4 

L I Ul ! TEST ) l J P 1 { 9 5 ) SLJI96) 

95 WRITE! 19, 133) 

133 FORMAT!* EXECUTION TERMINATED* /) 

RETURN 

96 UNLOAD 3 
ENDFILE N5 
ENDFILE N5 

RETURN 

97 WR I TE ! 9 , 100 l ) 

1001 FORMAT!* ERROR IN PROGRAM INPUTS* END OF JOB*} 
l CRASH - -1 
N5 = 9 
GO TO 46 

7CST SL J ( » ) ENIO(O) 

♦ THS(KBAD) SL J{ 39 ) 

NL 0 ! J ) = NLO! J ) + 1 
LO A ! KBAQ ) SLJ140) 

39 + THS ! NB AO } SL J ( S2 ) 
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SL J ( 7CS T ) 

40 ST A3 ( DATA ) 

SL J ( 7C ST ) 

92 NHIIJI * NHKJ) ♦ 1 

LDAINBAU) SL J I 40 ) 

300 SL J I * ) .UNPACK 'TIME* CHANNEL INTO X 

305 CALL READANY13, 1» DATAdU DATA! NWO ) * NWJ 
I F l NW ) 98, 306, 307 

98 LDAIN90) SAIML+1) 

SLJ (*) 

306 UNLOAD 3 
PAUSE 1 2 345 
GO TC 305 

307 I = [START 
L = LST ART 

CO 310 K=l,NB 

L0A4 (SHIFT) $ AL ( 3 10 A } 

31 0 A LDC1 ( DATA ) LLS(*J 

SCLI 7777777777777000B ) S T A3 ( X J 

1 = 1+ GELID 
310 L ^ POSID 

GO TC 300 

320 SL J t * ) .SET UP UNPACKING, AS IN GRAYMAP 

ISTART = ( NCHAN + 4)/5 
LSTART = NCHAN ♦ 5*11 - ISTART) 

CO 325 L = I, 5 

NEXTPOS * NCHAN +■ L 
DELIL) = I NEXTPOS - 1 )/5 
325 POSIL) = NEXTPOS - 5*CEL(L> 

NWD = ( NCHAN »NS S + 4>/5 
CO TC 320 

2STCL SLJi*> LDAI2CLN). 

♦SSKIIVC+2) LDAI2CLR). CHECK FIR REVERSE POLARITY DATA 
STAI2CL) SLJI2STCL ) 

2CLN LDA6 1 DAT A ) SUB2IIVC). CLAMP NORMAL POLARITY DATA 

2CLR LDA2IIVC) SUB6 ( DA T A I . CLAMP REVERSE POLARITY DATA 

502 SSK ( IKEY+3 ) SLJI507). JUMP IF DESKEWING ONLY 

SLJI 505) 

507 LDA ( J SKEW ) AJP3I97). CONTINUE IF DESKEWING ONLY 

N = NST 
K =0 

DC 535 IS=NA,NQ1,KP 

DC 534 J=1,NC 
K-K + l 

N = N +■ ISKEWI J) 

LDA6ICATA) STA3I DATA ) 

534 CONTINUE 

535 N - N + NKP 
GC TO 505 

END 
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PROGRAM AUTOCAl 

C USES PROCESS. NEW FSKIP, COMMENT, UNPACK 3 
C INPUT ON 3, OUTPUT ON A 

C CORRELATION REFERENCE CHANNEL IS A TYPED INPUT 

C ANSWER ‘LINE NO. =• WITH NS A NSB K5 NA MB BY 515 

C ANSWER 'CORRELATION REGION =' WITH Cl C2 BY 215 

C THE CORRELATION REGION SHOULD HAVE ELBOW ROOM INSIDE NA NB 
C TC PERMIT CORRELATION, AND NA NB SHOULD HAVE ELBOW ROOM 
C INSIDE 1 NS S TO PERMIT SLEW CORRECTION. 

C A PREVIOUS GRAYMAP OF THE CALIBRATION DATA IS VERY DESIRABLE. 

C SENSE SWITCH 2 ON TO PRINT CHANNEL CORRELATION INFORMATION 
C OUTPUT CF THE PROGRAM IS SLEW-CORRECTED AND DARK-LEVEL CORRECTED 
C AVERAGE LINE, DARK LEVEL, AND INTEGER AND FRACTIONAL CHANNEL SKEW 
C VIDEO DATA USED FOR SLEW COMPUTATION IS NOT DARK LEVEL CORRECTED 
C TO TURN OFF SLEW CORRECTION, MAKE Cl .LT, 0 
C IF KS IS -l, EVERY FIFTH ICHAN POINT ON N SA WILL BE TYPED 

DIMENSION DEL I 16 > ,N0FF( 16), OFFSET! 16 ) . YI400), PROD ( 1 01 ) , Z 1 2 000 ) , 
1CARKI 16 I 

COMMON DATA I 6000 I , X12000) 

COMMON RESERVE! B ) , NF, NR, MR, L0O, ID1, ID2, BANG, DANG, L90 
COMMON KEY, NPTS, NLINES, NSA, NSB, KS, NA, NB, KP, IS, TITLEI12), 
1 TAG1I12I, T AG 2 1 12 J , NSS, NCHAN, IDA, CONV, IPOS, IPACK, I NT 
EQUIVALENCES, Z), (Y.PROC), I Y I 102 I ,NOFF ! , <Y(115), OFFSET! 

INTEGER X, Y, DATA, DEL, DELO, DELI, DEL 2 , DEL3, TEST, FIRST, AVE, 
1 AVEMAX, SUM, SUMO, Cl, C2, CX, C2X, CD, DAT A I , PROD, PRODMAX 
SET INOEX K 1 =L» K2 = M» K 3 = N 
CAT A I NSKEW = 10) 

1 WRITE!9, 101 I 

101 FORMAT t*CCRRELAT ION REFERENCE CHANNEL =* > 

RE AO 19, 102) ICHAN 

102 FORMAT (815) 

WR I TE ( 9 , 106 I 

106 FORMAT ( *ARE YCU USING CALIBRATION DATA...*) 

READ (9, 107) ICAL 

107 FORMAT ( A8 ( 

IFIICAL.EQ.3HYESI 6, 8 

6 ICAL = 1 

GO TO 21 
8 ICAL = 0 

C READ AND PRINT TITLE AND OTHER INFORMATION 
21 CALL PROCESS! DATA, 1, 1) 

C READ ANO PRINT LINE AND COLUMN NUMBERS 

2 CALL PROCES S ! DA T A , 2, 1) 

AVAL = NPTS»NCHAN 
IFINVAL.GT.2000) 4, 5 

4 WRITE(9,104) 

104 FORMAT ( 20HNPTS*NCHAN .GT. 2000) 

GO TO 2 

5 CALL ZERO I X , X(NVAL) ) 

NY = NB - NA * 1 



I F (NY • GT . 40 0 ) 

58, 59 


58 

WR I TE ( 9 * LS8) 



158 

FORMAT ( *N6 - 

NA ♦ 1 

• G T « 


GO TO 2 
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59 WRITEI9.103) 

103 FORMAT! ‘CORRELATION REGION = *) 

RE AD ( 9 , 102 ) IC1, C2 



Cl = 

1 ABS ( I C 1 > 


I F < C 1 

.EQ.O) 12, 13 

12 

Cl 

= M I NO ( M A + 3 , NB-3J 


Cl 

= MAXOlClt NA) 


C2 

= Nfl - 3 

13 

C2 = 

M A X 0 < C 2 , Cl+1) 


C2 = MIN0IC2, NB ) 

V)RITEI4,167) Cl, C2 

167 FORMAT 1 6X ‘CORRELATION REGION =* 17, 15/) 

CX = ICl-NA I ‘NCHAN «■ ICHAN 
AC * C2 - Cl 

C2X = <C2-NA)*NCHAN ‘ ICHAN - NCH AN 
MX = MINOICl-NA, NSKEW) + 1 

KINC = NCHAN 

KS TART = I N A~1 ) ‘NCHAN 

NA2 « MAXOI NA, 3) 

NB2 = M I NO I NB , NSS-21 
FIRST = 1 
I NT = l 

NRANGE = MINOINSKEW, Cl-NA) + MINOINSKEH, NB-C2 ) + 1 
WK1TE<4,105> 

105 FORMAT ( *0 LINE NA NB*/) 

C UNPACK A LINE OF INTEGER DATA INTO THE ARRAY DATA 
3 CALL PROCESS ( DAT A , 3) 

IFIKS.LT. 0) 400 , 71 
71 I F ( KEY . EC .4 ) 35, 75 

C ESTIMATE AD JAC ENT- PO I N T DIFFERENCE ON THE HIGH SIDE 
75 LOA(FIRST) AJP(U) 

CO 10 J * 1, NCHAN 
K 1 = KS TART + J 
K2 = K1 + KINC 
CO 7 I = 1, NPTS-1 

Y(I) * IABS(DATA(K2)-DATA(K1) ) 

K1 - K l + KINC 
7 K2 = K2 + KINC 

CALL SORT 1 1 Y, NY, -1) 

10 CELU) = Y(4*NPTS/5) * 2 

C EDIT NOISE SPIKES CUT OF THE CATA 

11 CO 20 J = 1 , NCHAN 

K 1 = ( NA2-1 ) ‘NCHAN + J 

KO = K1 - KINC 

K2 = Ki + KINC 

K3 - K2 + KINC 

CELO = DATA ( KO) - DAT A( KO-K INC ) 

CE LI = DAT A ( K 1 ) - DATA(KO) 

CEL2 = DAT A I K2 ) - DATA(Kl) 

CO 25 I = NA2, NB2, KP 

DEL 3 = DATAIK3) - DATAIK2) 

LDAIDEL1I SCM ( DEL 2 ) AJP2I24) .IF DELI, DEL2 SAME SIGN, TCI 24 
TEST = M I NO I I ABS I DEL 1 ) , IABSIDEL2)) 

LEVEL = MAXOI IABSI DELO) , l ABS ( DEL 3 ) , DELIJU 
IF I TEST .GT. LEV EL *5) 22, 24 
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22 K1 = K2 - KINC 

DAT A ( K 1 1 = ( DATA! Kl+NCHAN ) + DAT A ( K 1-NCH AN )) /2 

24 DELO = DELI 

DELI * DEL2 
DEL2 = DEL 3 

K2 = K2 + KINC 

25 K3 = K3 + KINC 

20 CONTINUE 

C CORRECT FOR SLEW OF THE CALIBRATION DATA 

LDA(ICI) A JP 3 ( 6 2 3 
LOAUCAL) A J P { 62 ) 

LDA(FIRST) A JP 1 ( 31 3 

IF INL INES.EQ.2 .OR. NLINES.EQ.6) 60, 65 

60 L = CX 

NL = NLINES - 1 
CO 61 K = 1, NC 
M = L + NCH AN 
V ( K ) = |X(M-xa)l/NL 

61 L - M 

65 CD = KST ART ♦ CX 

SLJM200). SUP Y ( K ) * DEL DAT AIL) FOR ICHAN AT CAL REGION 
SUMO = SUM 
CD = CD ♦ NCH AN 
SLJ41200) . 

IF (SUM.GT .SUMO ) 26, 27 

26 SUMO = SUM 

NA = NA + 1 

NB = NB + 1 

IF(NB.GT.NSS) 56, 28 

56 WRITE(9,157> NA, NB, NSS, IS 

157 FORMAT (^CALIBRATION HAS SLEWED OFF END. NA =* 14, *, NB =* 14, 

l *, NSS =* 14, *, LINE NO.* 153 
GO TO 2 

28 CO * CD ♦ NCHAN 
SL J4 ( 200 ) . 

IFISUM.GT.SUM0> 26, 31 

27 CD = CD - NCHAN 

33 CD = CD - NCHAN 

S L J4 ( 200 J . 

IF(SUM.GT.SUMO) 32, 31 
32 SUMO = SUM 

NA = NA - 1 
NB = NB - 1 
IFINA.LT. 1) 56, 33 
31 K ST ART = ( N A- 1 ) «NCHAN 
NA2 - MAXOlNA, 3) 

NB2 = MINOINB, NSS-2 J 
IFISENSE SWITCH 2) 69, 70 
69 WRITE(9,131 > NA , NB 

70 WR ITE ( 4 , 131 3 IS, NA, NB 
131 FORMAT (316) 

C UPDATE THE SUMS OF THE DESIGNATED POINTS 

62 FIRST = 0 
K = KSTART 

CO 29 L = 1, NVAL 
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K = K ♦ 1 

29 X ( L ) - X ( L ) + DAT A ( K ) 

I F ( K E Y . 6Q. 2 ) 3, 35 

C AUTOMATIC DESKEWING CALCULATION 

35 LC = MAXO ( C X - ICHAN - NSK E W*NCHAN, 0) 

DO 50 J = 1, NCHAN 

IF< J.EQ. ICHAN) 50, 40 

40 PRODMAX = 74COQOOOOOOOOOOO0 
LC J - LC + J 

CO 48 M » 1 9 NRANGE 
PROD ( M ) = 0 
L * LC J 

DO 47 K - CX, C2X, NCHAN 
N = L + NCHAN 

PROO(M) = P ROD ( M ) ♦ ( XI K*NCHAN)~X<K )) * (X(N)-X(L)> 

47 L = N 

1 F < PROD I M ) • GT • PRODMAX) 41, 48 

41 PRODMAX - P ROD ( M ) 

M2 = M 

48 LC J = LC J ♦ NCHAN 

IFIM2.EQ.1 .CR. M2 , EQ . NRANG E ) 42, 43 

42 H I AS = 0. 

GO TO 44 

43 Ml = M2 — 1 
M3 = M2 + 1 

BIASN = PROD ( M3 J - PR0DIM1) 

BIASD = 2 * ( 2* PROD (M2) - PROD(Ml) - PR0D(M3>) 

BIAS = BIASN/BIASO 

44 NOFFIJ) * M2 - MX 

OFFSET (J) * FLO A T ( NOFF I J ) ) + BIAS 
I F ( $ EN5 E SWITCH 2J 30, 50 

30 WRITE(4,150) BIASN, BIASD, BIAS, CX, C2X, LC, LCJ, NRANGE, MX, 
1 Ml, M2, M3, PRODMAX, (PROD(M) f M" 1 , NRANGE ) 

150 FORMAT ( 1H0 3F12-4, 918/(0115)) 

50 CONTINUE 

NOFF ( ICHAN) - 0 
CFFS ET ( ICHAN) * 0. 

C ESTABLISH DARK LEVEL FOR EACH CHANNEL 
FN = NUNES 
FN = 1 ■ / ( FN *CONV ) 

LDA(ICAL) A JP ( 72 ) 

CO 55 J - It NCHAN 
SLJ41300). COMPUTE DARK LEVEL 

55 DARK ( J ) = FLOAT (AVEMAX)/8* * FN 

WRITE (4, 155 ) ( DARK ( J ) , J*l» NCHAN > 

155 FORMAT ( *0 DARK * 12F8.3) 

WRITEI4, 156) 

156 FORMAT (1H > 

C COMPUTE AVERAGE POINTS SUBTRACTED FROM THE DARK LEVEL 
K - 0 

DO 36 l = 1, NPTS 

DO 36 J = 1, NCHAN 

K = K + 1 

36 Z < K > = DARK(J) - FLO AT i X ( K ) ) *FN 
GO TO 74 
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72 DO 73 K * 1 ( NVAL 

73 ZIKI - FLOAT I X ( K ) I »FN 

74 K 1 = 1 

K2 = NCHAN 

CO 37 I = NA, N8, KP 

WRITE(4,137) I, I Z I K > i K=Kl f K2l 
137 FORMAT I 16 » 12F8.3) 

K1 = K1 ♦ NCHAN 
37 K2 = K2 + NCHAN 

WRITEI4.144) INCFF(J), J=l, NCHAN) 

144 FORMAT ( *OSKEW« 15, 11181 

WR IT E ( 4 , 145 I (OFFSET ( J ) , J=l, NCHAN) 

145 FORMAT [ 6X 12F8.3) 

GO TC 2 

200 SLJ<») .SUM CF Y I K I * DEL DAT AIL) FOR [CHAN AT CAL REGION 
SUM = 0 

L = CD 

CO 201 K = 1, NC 
M = L ♦ NCHAN 

SUM * SUM ♦ Y(K ) *( DATA! M )-DATA< L I ) 

201 L = M 
GO TG 200 

300 SL J ( * ) . DARK LEVEL COMPUTATION 

AVE = 0 
L = J 

CO 51 I = 1, 8 
AVE = AVE + X(L) 

51 L = L ♦ NCHAN 

AVEMAX = AVE 
K = J 

CO 54 I * 9, NPTS 

AVE = AVE - X ( K ) + XI L) 

K * K + NCHAN 
L = L + NCHAN 

54 AVEMAX = MAXO ( AVEMAX , AVE) 

GO TC 300 

C PROGRAM TO TYPE EVERY FIFTH ICHAN VALUE ON LINE NSA 
400 L = ICHAN 

LINC = 5 *NCH A N 
LINC2 = LINC + LINC 
DO 401 I » 1, NSS, 10 

WRITE(9,402) I, DATA(L), DATA ( L«L INC ) 

402 FORMAT (13, 215) 

401 L = L + LINC2 

GO TC 2 
END 


PROGRAM EXPMAP 
C VERSION 1.0 2/8/72, ZUK 

C NEEDS MAPTRAN, 4-ARG BUFPACK, PROCESS AND UNPACK3 
C LIBRARY ON 1, CONTROL DATA ON 2, ADTESTED SCANNER DATA ON 3, 

C 2-CHANNEL RECOGNITION AND EXPONENT OUTPUT ON 4, HISTOGRAMS ON 5 
C PROVICES FOR MAPPING WITH PREPROCESSED DATA. FOR UNTRANSFORMEO 
C DATA, NTRAN = 1 AND DATUM! J) = CORRECT SUBSET OF DATA! I » J ) . 

C NT = NO. OF SIGNATURES 

C NO = NO. OF CHANNELS ON TAPE 

C NN = NO, OF CHANNELS IN SIGNATURE DECK 

C NV = NO. OF CHANNELS IN SIGNATURE USEO 
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C SO FT = NO. OF SQ FT PER RESOLUTION ELEMENT 
C S(I> = INTEGER REPRESENTING SIGNATURE I 

C TAG! 1 > « T AGA I I ) = 16-CHARACTER IDENTIFICATION FOR HISTOGRAMS 
C I SUB = SIGNATURE DECK CORRESPONDS TO A SUBSET OF DATA CHANNELS 
C CODEC J ) DEFINES THIS SUBSET 

C JSUH = USED SIGNATURE IS A SUBSET OF THE SIGNATURE DECK 
C ICODEIJI DEFINES THIS SUBSET (ORIGINAL CHANNEL NUMBERS) 

C NTRAN = NO. OF TRANSF'NS COMPUTED. UNTRANSFORMED COUNTS AS 1. 

C NTRAN NEGATIVE TO SIGNAL A SPECIAL OPERATION 

C I MENS I ON AMAPI 6000 ) > HI2280), ITABLEI52) 

COMMON AI20.12), Btl2.12.20), D(20l, DATAI6000), E IGt 12) . 10(10), 

1 IBINI21.52), ICI12), ICOCEI 12) . IPT(IOO), NPt21), S(2l), S2I21), 

2 T AG ( 21 ) , T AGA (21), Z(12> 

COMMON CODE 112), 0ATUMI12), NN, NTRAN 

COMMON R(S), BAN2, 0AN2, CC» NF, NR, MR, L80, ID1, ID2, BANG, DANG 
COMMON L90, KEY, NPTS, NLINES, NSA, NSB, KS, NA, NB, KP, IS, 

1 TITLEI12), TAG1I12), TAG2I12), NSS, NCHAN , KR, CONV, IPOS, 

2 IPACK, INT 

EQUIVALENCE IB, HI, (DATA, AMAP ) 

INTEGER CODE 

INTEGER AMAP, REPLY, S, TAG, TAGA, TAG1, TAG2, TITLE 
INTEGER SEC, SEC60 

CATA I IT ABL E( 11= 511, 490, 480, 470, 460, 450, 440, 430, 420, 410, 


1400, 

390. 

380, 

370. 

360, 

350, 340, 330, 

320, 

310, 300, 290, 

280 

2270, 

260. 

250, 

240 , 

230, 

220, 210, 200, 

190, 

180, 170, 160, 

150 

3140, 

130, 

120, 

110, 

100, 

90, 80, 70, 60, 

50, 

40, 30, 20, 10 

, 0) 


LIB(BUFWAIT) 

LOCI TEST * 70) 

PARTMAP 

20 FORMAT (2015) 

21 FORMAT! 10A8) 

22 FORMAT I 5E15 .8 ) 

SLJH504) 

CALL CORECON 
CUMMY = 0. 

501 CALL PROCESSIDATA, 1, 0, 3, 2) 

QRITEI5, 1000) TITLE 
CONV 10 = CONV/IO. 

503 CALL PROCESS ( DATA, 2) 

ENA< BUFNA IT) SALIL80) 

CALL 8UFF0N 
IF IMR.GE.O) 17, 642 

17 RE AO ( 2, 200) NT, NO, NN, NV, LTRAN, SQ FT 

200 FORMAT(5I3,3X,5E12.0) 

EXPLIM = 99.6 

622 FKS = KS 
FKP = KP 

ACRE = SC FT * FKS * FKP / 43560. 

NT1 * NT ♦ 1 
I SUB = NO - NN 
JSUB = NN - NV 
NTRAN = IABS(LTRAN) 

IF 1 1 SUB ) 305, 306, 305 

C * READ A SUBSET OF DATA CHANNELS IF DESIRED. 

305 RE ADI 2 , 20 ) ICODE(J), J = l,NN) 

CALL SORT 1 1 CODE , NN, -1) 

GO TO 308 
CO 307 J = 1 , NN 


306 
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307 CODE < J ) = J 

300 1 F ( J SUB *GT * 0 ) 330, 331 

C * READ A SUBSET OF SIGNATURE CHANNELS* 

330 RE AD < 2 , 20) IICCDEU), J=1,NV) 

CALL SORT 1 < I CODE , NV * -1) 

GO TO 337 

331 DO 332 J = 1* NV 

332 I CODE ( J ) = CODE ( J ) 

337 L = 1 

CO 340 J * I* NN 

IF( ICODECL) .EQ.CODEU)) 341, 340 
341 I C { L ) = J 

L = L + 1 
340 CONTINUE 

314 S(NT*l)=5ll 

TAG< NT+ 1 ) = 8HNCT CLAS 
TAGAtNT + 1) = 6HS I FED 

7 L S T ART = 1 

C * * # * READ SIGNATURE DECKS AND ASSOCIATED INFORMATION * * * * 

CO 108 I - 1 , NT 
READ < 2 » 21 ) ID 

IF( ICU) *EQ.4HSAME) 108, 55 
C READ SIGNATURE CECK 

55 REACC2, 22) ( ( A (I , J > , J= 1, NN ) , { ( B ( J , K , I > , K- J , NN ) , J= 1 , NN ) ) 

C READ INTEGER VALUE ASSOCIATED WITH SIGNATURE, NAME OF SIGNATURE, 

RE AD 1 2 , 600) SCI), TAG(I), TAGA(I) 

600 FORMAT C I5,2A8, 3F12.0) 

C REARRANGE THE MATRIX ACCORDING TO THE SUBSET* 

104 CQ 2 J= 1 , NV 

JJ = ICOCECJ) 

CO 2 K=J,NV 
KK = I CODE ( K ) 

B ( J , K , I ) = BUJ,KK, I > 

2 B ( K , J , I ) = 61 J,K, I ) 

C CHECK TG SEE IF THE MATRIX IS POSITIVE DEFINITE. 

IFCDIAG(B(OtltI), EIG, DUMMY, -NV, 12)) 23, 334, 23 
23 WRITE(9,124) TAG(I), TAGA(I), (EIG(J), JM,NN) 

124 FORMAT ( *M ATR I X * 2A8, • IS NOT POSITIVE DEFINITE. EIGENVALUES * 

1 *ARE */ ( X 2£ 10 * 3 ) ) 

STOP 

334 CALL INVERSEiBi 0,1, I ), 6(0, 1,1), NV, 12, 12, DU), ISERR) 

WRITEI5, 123) TAGII), T AG A { I ) , SCI), DC I > 

123 FORMAT ( 2A8 , 15, E15.5) 

336 C ( I ) = A LOG ( G I I )) 

C 

C ARRANGE THE MATRIX B SO THAT IT FITS CORRECTLY INTO THE LINEAR 
C MATRIX H* 

L = LST ART ♦ 1 
H ( LST ART ) = B ( 1 , 1 , 1 ) 

DO 105 J=2, NV 
Jl = J - 1 
CO 107 K= 1,JX 

HU) = B( K, J, I ) « 2. 

107 L s L ♦ 1 

H ( L ) = 6 ( J , J , I ) 


70 



2eri 


FORMERLY WILLOW RUN LABORATORIES. THE UNIVERSITY OF MICHIGAN 


105 L = L + l 

108 LSTART = LS T ART + 144 

CALL ZERCINPI 1 ) ,NP( NT + 1 )) 

C 

C ** OUTPUT CATA IS ASSUMED TO BE PACKED AND POSITIVE. 

NCHA1 a 2 
NSS1=!NB-NA l/KP + l 

TAGl(l) = BHRECOGN I T 
TAG2I1) = 8HICN MAP 
T AG1 ! 2 ) = 8HEXP0NENT 
TAG2I2) « 8H / 10 

HRITE14) TITLE, NSSl, BAN2, DAN2 » NCHA1 , CONV, 1, 1, TAG1, T AG2 

652 CALL ZERC! 1 BIN, IB1N! 1092) ) 

CALL TIMER! 3HSET ) 

C READ AND UNPACK THE INPUT DATA. 

642 CALL PROCESS! DATA, 3) 

LDAIKEY) INAI-4) AJP ( 51 A ) 

CALL TIMER! 2HG0 I 

1 = 1 
N = Q 

W0R0=0 

DO 10 IP = NA, NB, KP 
IPP=< I P-1) *ND 

C MAPTRAN TAKES THE CATA I AND IF NECESSARY TRANSFORMS IT) POINT BY 
C POINT AND PUTS IT INTO THE SMALL ARRAY DATUM. 

CALL MAP TRAN! DATA I IPP > ) 

N=N+ 1 

LST ART * 1 
GM IN = 1.E10 
IMIN = NTl 

C * « RECOGNITION CALCULATIONS * » 

DO 12 I M= l , NT 
L = LST ART 
JJ = IC11) 

Z ( 1 ) = DATUM I J J I - AIIM.JJ) 

SUM = Z(l) • 2111 * HCL ) 

L = L + 1 
C 

CO 11 J=2,NV 
JJ = 1C!JI 

Z ( J ) = DATUM! JJ) - A I IM, J J ) 

SUM2 = 0. 

DO 13 K = 1 , J 

SUM2 = SUM2 ♦ Z!K)*HIL) 

13 L = L ♦ 1 

11 SUM = SUM + Z ! J ) *SUM2 

IFISUM.GT.EXPLIM) 12, 24 
24 GM = SUM ♦ D( IM ) 

IF(GM.LT.GMIN) 14, 12 

14 1 M IN = IM 
GMIN = GM 
GSUM = SUM 

12 LSTART = LST ART ♦ 144 
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C AM A P IS THE OUTPUT ARRAY 
AM AP ( N ) = S I I MIN) 

N = N* 1 

AMAPIN) = G SUM*C0NV 1 0 + .5 

C SEARCH FOR MAKING A HISTOGRAM. 

EN12<52) 

♦ THS2 ( I T ABL E ) SLJ1900) 

900 J = J+1 

IBINUMIN.J) = I8I.N! IMIN.J) + 1 
C NP * NUMBER OF POINTS RECOGNIZED FOR EACH SIGNATURE. 

NP(IMIN) * NPIIMIN) + 1 
10 CONTINUE 

C WRITE THE DATA. 

CALL BUFPACK ( AMAP, N, AMAP, 4) 

CALL TIMER! 4HSTCP) 

C LINES NOT FINISHED GO TO 642 
C LINES FINISHED CO TO 51 

IFIKEY .EQ. 2) 642, 51 

51 A NS6 = IS - 1 

51 WR1TEI9, 1022) 

1022 FORMAT!* CONTINUE FILE...*) 

RE AD ! 2 , 1021 ) REPLY 
1021 FORMAT! Al, 12A8) 

IF ( REPLY .EG. 1HY) 503, 504 
C •* YES OR Y GO TC 503 
504 ENDFILE 4 

CALL BUFF OF F 

CALL TIMER! SEC, SEC60 ) 

IHRS = SEC/3600 

MINTS = (SEC - 3600*1 HR S)/60 

I SEC = (SEC -3600* IHRS - 60*MINTS) 

NL1NES = ( NSB - NSAI/KS + 1 

WRITEI5, 780) NLINES, NSS, ND, NPTS, NV, NT, IHRS, MINTS, I SEC 
780 FORMAT!/* INPUT...* 110* LINES, *15* P TS/L I NE , * I 5* CHANNELS*/ 

1 * PROCESSED. ..*15* POINTS, *15* CHANNELS ,* I 5* SIGNATURES*/ 

2 10X *TCT A L TIME* * 15* HRS,* 15* MINUTES, *15* SECONDS*//) 

80 WR I T E ! 5 * 32) 

32 F0RMAH//7X, 4 HNAM£,14X,* VOLTAGE NUMBER ACREAGE*/) 

DO 72 J=1 , NT+ 1 
FNP*NP ( J ) 

AREA=FNP*ACRE 
FS=SI J) 

S2 ( J I =FS/C0NV 

WRITE ! 5 , 81) T A G ! J ! , TAGA(J), S2!J)» NP(J), AREA 
81 FORMAT! 5X 2 A8 , 3X , F4. 1 , 1 1 2, F 12. 2 1 

72 CONTINUE 

C * » » HISTOGRAM * * * 

651 DO 660 1=1, NT 
FNP= NP ( I ) 

AREA=FNP«ACRE 

C WRITE HEADER 

WRITEI5, 1000) TITLE 
1000 FCRMATllHl, 10X,12A8) 
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WR I T £ { 5 , 1001) T AG ( I ) , TAGA(I), S2U), NPU)» AREA, DU) 

L001 FORMAT ! / / 7X 4HNAME 14X •VOLTAGE NUMBER ACREAGE LN(DET)*/ 

1 5X 2 A3 , 3X FA * 1 112, 2F12.2) 

FT AB = ( F LC AT ( ITABLEl 1) ) ) /CONV 10 
GTAB = (FLOAT ( ITABLEI51) ) ) / CON V 1 0 

WRITE<5, 1002) FTAB, lBtN(I,l), GTAB, IBINU,52i 
1002 FORMAT ( 1 H *NC. CF POINTS .GT.*F5.1* =* 18 *, NO* OF POINTS .LT.* 
1 F5 * 1 • = * 18 //* EXPONENT* 2H**) 
l JK = 0 

C DETERMINE THE BIN WITH THE LARGEST NUMBER OF POINTS 
DO 662 J=2, 5 l 

662 UK = MAXOIIJK, IBIN(I,JJ) 

EN 1 2 ( 50 ) 

C CALCULATE PERCENT IN EACH BIN. 

664 IJP21665) SLJ(660) 

665 PCT = (FLOATI I B I N ( I , J+2 ) )/FL0AT(NP{ I ) ) )*100. 

K= ( 1 00* I B IN ( I , J +2 ) ) / I JK 

CALL ZERO! IPT( 1 ), IPT( 100), 1H ) 

IF(K.GT.O) 66 T, 66B 

C FILL ARRAY AND WRITE* 

667 CALL ZERCUPTU), IPT(K), 1H* > 

668 FTAB = I FLOAT U TABLE ( J+2) M /CQNV10 
GTAB * (FLOAT! ITABLE(J+1) ) )/C0NV10 
WR I TE ( 5 , 1005) FTAB, GTAB, < 1PT IK) , K«1 , 100 ) » PCT 
1005 FORMAT ( 1H F4.1 * TO * F5*l, 100A1 F6.2) 

GO TO 664 

660 WR I TE { 5 i 2006) 

2006 FORMAT!// 4H *•* CHANNEL 2 IS EXPONENT / 1 0* * ) 

65 1 F ( REPLY • E Q. 1HD) 506, 501 

C IF CONE. GO TO 506 . . *0T HER W I SE , GO TO 501 

506 ENDFILE 4 
ENDFILE 4 

LIUl(TEST) IJPK721) 

ENDFILE 5 
ENDFILE 5 
UNLOAD 2 
UNLOAD 3 
UNLOAD 4 
UNLOAD 5 

721 WRITE! 19,1030) 

1030 FORMAT ( • EXECUTION TERMINATED* ) 

STOP 
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